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ABSTRACT 


This  report  summarizes  the  work  completed  by  May  31,  1981  for  the  Air 
Force  Office  of  Scientific  Research  as  part  of  Grant  No.  AFOSR77-3336.  The 
research  deals  with  the  broad  topics  of  initiation,  combustion  and  transition 
to  detonation  in  homogeneous  and  heterogeneous  reactive  mixtures.  One  work 
area  is  involved  with  certain  aspects  of  ignition  source  effects  in 
reactive  fuel-air  mixtures*  These  aspects  include  effects  of  chemical 
sensitizers,  flame  acceleration,  flame  area,  and  ignition  point  location. 

The  other  area  involves  the  hydrodynamic  modeling  of  ignition  and  flame¬ 
spreading  in  granular  energetic  solids  to  predict  the  potential  for 
deflagration-to-detonation  transition  (DDT). 

Key  results  in  the  first  area  are  that  chemically  sensitized  clouds 
can  lead  to  detonation,  that  flame  acceleration  or  a  large  increase  in 
the  time  rate  of  increase  of  the  flame  area  are  needed  for  transition 
from  deflagration  to  detonation  and  that  it  is  very  difficult  to  generate 
damaging  overpressure  from  edge-ignited  combustion  even  for  very  high 
subsonic  burning  velocities. 

The  research  dealing  with  analysis  of  DDT  in  porous  high  energy  solid 
propellant  has  shown  (for  the  first  time)  actual  steady-state  detonation 
solutions,  following  the  unsteady  flow,  for  materials  with  sufficient 
porosity  and  critical  burning  rate  properties.  Limits  of  the  run-up 
length  to  detonation  are  predicted  as  a  function  of  propellant  chemical 
energy,  burning  rate,  bed  porosity,  and  granulation  (size).  The  detonation 
states  conform  to  realistic  measured  conditions  for  porous  HMX  and  RDX 
propellants. — 
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RESEARCH  OBJECTIVES/ STATEMENT  OF  WORK 


The  research  deals  with  the  broad  topics  of  initiation,  combustion 
and  transition  to  detonation  in  homogeneous  and  heterogeneous  reactive 
mixtures.  Specific  areas  of  research  include  the  following: 

A)  Analytical  work  directed  to  study  direct  initiation  of  detonation 
by  a  nonideal  blast;  wave  in  chemically  sensitized  reactive  fuel-air  clouds. 

B)  An  experimental  study  of  non-linear  effects  upon  initiation  in  a 
reactive  fuel-air  mixture  with  initiation  produced  behind  the  Mach  stem 
shock  generated  on  a  ramp. 

C)  An  experimental  study  of  the  effect  of  obstacles  on  flame  acceleration 
in  cylindrical  geometry  and  in  layered  clouds. 

D)  The  numerical  and  theoretical  calculation  of  the  blast  waves 
produced  by  a  centrally  ignited  spherical  source  region  and  application  of 
acoustic  monooole  source  theory  and  numerical  techniques  to  the  deflagrative 
and  detonation  combustion  of  a  cloud  of  arbitrary  shape. 

E)  The  DDT  phenomenon  in  granulated  propellant  or  explosives  involves 
a  series  of  complex  transient  processes  that  are  not  well  understood  at 

the  present  time.  It  is  hypothesized  that  the  normal  burning  process  of 
the  solid  propellant  is  disturbed  by  an  abnormality  such  as  a  crack  in  the 
propellant  grain.  This  abnormality  generates  regions  of  porous  propellant 
which  can  be  ignited  locally,  causing  a  pressure  buildup  and  formation 
of  a  weak  shock.  If  detonation  is  to  be  excited  following  this  ignition, 
it  is  necessary  to  ensure  a  sufficiently  rapid  pressure  buildup.  In  the 
case  of  porous  propellant,  this  may  be  achieved  as  a  result  of  the  penetration 
of  gaseous  combustion  products  into  the  interior  pores  of  the  solid,  which 
leads  to  the  disturbance  of  surface  burning  conditions.  Thus,  in  this  case, 
heat  transfer  by  conduction  is  replaced  by  convective  heat  transfer. 

Subsequent  acceleration  of  ignition  (flame)  fronts  begins  and  pressure 
waves  are  generated  which  become  shocks.  These  shocks  cause  large  local 
over-pressurization  and  often  change  into  a  detonation. 

To  analyze  this  phenomenon,  the  reactive  two-phase  (solid,  gas)  con¬ 
servation  equations  of  continuity,  momentum  and  energy  must  be  solved  along 
with  any  constitutive  relations  to  account  for  heat,  transfer  interaction, 
pressure  losses  through  the  aggregate,  ignition  criteria,  unsteady  burning 
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rates,  etc.  It  is  only  though  solutions  of  such  a  fluid-mechanics  model 
that  one  can  develop  criteria  to  specify  the  conditions  under  which  the 
burning  propellant  is  susceptible  to  transition  to  detonation.  Considerable 
effort  is  involved  in  developing  advanced  numerical  integration  schemes  to 
adequately  handle  the  prediction  of  strong  shock  waves  in  such  transient 

flows . 

From  the  solution  of  the  unsteady  two-phase  flow  model  it  is  possible 
to  assess  the  important  physical  and  chemical  factors  that  are  inherit 
in  the  DDT  phenomenon.  It  is  then  also  possible  to  interpret  the  experiments 
where  defiagration-to-detcmation  transition  occurs,  to  assess  whether 
such  factors  as  granulation,  compaction,  confinement,  chemical  content, 

and  burning  rates  can  be  properly  tested. 

Experimental  research  to  evaluate  the  necessary  constitutive  relations 

for  the  momentum  and  heat  transfer  in  porous  beds  at  high  flow-Reynolds 
number  :.s  also  underway,  since  these  data  represent  input  to  the  modeling 

effort. 
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STATUS  OF  THE  RESEARCH  EFFORT 

A.  Direct  Initiation  of  Detonation  with  Chemically  Sensitized  Clouds. 

The  question  to  be  answered  by  numerical  calculation  was  whether 
the  direct  initiation  of  detonation  could  be  generated  by  an  initiating 
shock  wave  produced  by  a  proper  distribution  of  a  chemical  accelerator  in 
a  source  region  within  a  fuel-air  mixture.  The  reaction  was  modeled  by 
a  modified  Arrhenius  law  of  the  form, 


dA 

37 


[‘♦*2 


exp  (-E*/e) 


where  A  is  reaction  coordinate,  A1  is  the  dimensionless  prefactor  for  the 
reaction,  ^  accelerator  effectiveness  factor,  E*  is  the  dimensionless 

activation  energy,  v  and  e  are  dimensionless  specific  volume  and  internal 
energy  respectively  and  n  is  the  order  of  the  reaction.  The  function 
F (6)  is  given  by 


R_  -  R 


F(6)  =  cos  3u(5  -  9  cos  irfi  +  8  J/16 


and  R  is  the  radius  of  the  core  where  the  chemical 


where  6  *  = - = 

Ro-*l 

accelerator  is  uniformly  distributed  and  R^  is  the  outer  edge  of  the 
transition  zone  for  the  distribution  of  the  chemical  accelerator. 

A  modified  Lagrangian  time  dependent  finite  difference  (Qppenheim 
CLOUD)  program  was  used  in  the  study.  Early  efforts  using  the  reaction 


rate  law  Q  =  Q  w^ere  Q  is  the  total  energy  of  the  reaction  led  to  such 

*  * 

large  values  of  Q  that  it  was  necessary  to  limit  the  maximum  value  of  Q  in 

order  to  keep  a  sufficiently  large  time  step  so  the  computing  time  would 
not  become  excessive. 

Limiting  values  of  the  controlling  parameters  in  which  transition  to 
detonation  occurred  were  found  in  1973  for  this  model.  Subsequently  in 
1979  we  sought  to  investigate  the  effect  of  a  more  rapid  cutoff  in  reaction 
rate  below  1000  K.  Our  curve  for  the  logarithm  of  reaction  rate  versus 
inverse  temperature  puts  the  "knee"  at  1000  K  and  contains  a  parameter 
which  shifts  the  value  of  the  rate  at  the  knee.  As  we  reduced  the  value 
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of  this  parameter  tho  reaction  2one  became  so  narrow  that  it  encompassed 
only  one  cell.  Therefore,  we  sought  to  modify  the  CLOUD  program  by  sub¬ 
dividing  a  band  of  cells  which  moved  with  the  reaction  front.  A  linear 
curve  fit  led  to  numerical  calculation  instabilities  and  a  cubic  spline 
function  fit  permitted  some  extension  of  the  range,  but  we  have  not  been 
able  to  solve  the  problem  satisfactorily. 

Besides  the  annual  reports  to  the  contractor,  this  work  led  to  tho 
publication:  H.  0.  Barthel,  and  R.  A.  Strehlow,  "Diroct  Detonation 
Initiation  by  Local ized  Enhanced  Reactivity,1'  AIAA  Paper ,  79  0286  (1979). 

B.  Nonlinear  Effects 

The  nonlinear  initiation  experimental  program  in  1977-1978  used, 
as  a  geometry  for  reflected  shock  initiation  under  non-uniform  conditions, 
a  15°  ramp  at  the  end  of  the  shock  tube  placed  just  before  the  shock  hits 
the  back  wall.  The  Mach  stem  shock  would  generate,  when  reflected,  a 
temperature  behind  it  approximately  100  K  higher  than  the  bulk  gas  reflected 
temperature  due  to  the  normal  shock  reflection  of  the  incident  shock.  The 
initial  shock  strength  was  adjusted  so  that  the  normal  reflected  shock 
temperature  lies  in  the  range  from  approximately  900-1200  K  on  different 
shots,  and  ♦’erefore  the  ratio  of  time  delays  in  the  regions  behind  the 
reflected  normal  and  reflected  Mach  stem  shocks  could  be  adjusted  over 
a  large  range  because  of  the  nonlinear  behavior  of  the  slope  of  the  induction 
delay  curve  versus  1/T  in  this  temperature  regime.  The  program  was  less 
than  successful,  for  we  did  not  obtain  meaningful  smoked  foil  records  and 
had  to  resort  to  pressure  records  to  detei'mine  that  transverse  initiation 
occurred  on  some  runs.  The  results  appear  in  the  report  by  John  K.  Soldner, 
"Direct  Initiation  of  Detonation  Using  Finite  Amplitude  Wave  Acceleration," 
MS  thesis,  University  of  Illinois,  1978. 

C .  Propagation  of  a  Flame  in  a  Layered  Propane-Air  Mixture 

In  connection  with  the  development  of  monopole  source  theory  for 
the  deflagrative  combustion  of  a  cloud  of  arbitrary  shape  (this  will  be 
discussed  in  the  next  section) ,  it  was  necessary  to  determine  how 
rapidly  a  flame  propagated  in  a  layered  cloud  under  conditions  where  it 
was  edge- ignited.  For  this  purpose  a  80  cm  by  80  cm  aluminum  plate 
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was  mounted  at  a  10  degree  angle  to  the  horizontal  and  a  8  cm  by  8  cm  by 
20  cm  high  chamber  was  filled  with  a  combustible  propane  air  mixture  and 
mounted  on  the  top  e'ge  of  this  plate.  When  this  chamber  was  tilted  the 
propane-air  mixture  eoured  out  over  the  plate.  A  torch  was  mounted  at 
the  lower  edge  of  the  plate  and  when  the  propane-air  mixture  reached  the 
torch,  it  flashed  back  to  the  source.  Schlierer.  photographs  were  made  to 
determine  the  flame  propagation  velocity.  It  was  found  that  the  flame 
propagated  at  about  3  to  3^  times  the  maximum  normal  burning  velocity 
for  a  propane-air  mixture.  This  information  was  used  to  construct  the 
theory  which  will  be  discussed  in  a  later  section  (see  section  D) . 

Propagation  of  a  flame  in  a  2-dimens.ional  layered  mixture 

Preliminary  to  the  study  of  the  effect  of  obstacles  on  flame  acceleration 
an  apparatus  was  built  which  was  10  cm  wide,  50  cm  long  and  IS  cm  high 
containing  a  bottom  plate  of  aluminum  with  a  bronze  insert  to  allow  the 
infusion  of  combustible  mixtures,  two-side  glass  walls  and  an  endplate 
containing  5  spark  gaps;  for  ignition.  On  the  opposite  side  of  the  endplate 
a  gate  was  placed  which  was  used  to  control  the  layer  height  and  was 
dropped  out  of  position  just  before  the  sparks  were  fired.  This  apparatus 
was  used  to  determine  flame  propagation  behavior  and  the  effect  of  the 
buoyancy  vtctor  on  flame  propagation  using  methane-air  and  propane-air 
mixtures.  (Methane-air  is  lighter  than  air  and  therefore  the  apparatus 
was  mounted  upside  down  for  the  methane  tests).  The  preliminary  results 
of  these  tests  as  well  as  the  tests  on  flame  propagation  over  the  inclined 
plate  are  reported  in  a  master's  thesis  by  P.  G.  Huseman,  dated  1980, 
entitled,  "Two  and  Three  Dimensional  Unconfined  Flame  Studies." 

D.  Theoretical  Analysis  of  the  Blast  Wave  Produced  by  Deflagrative  and 
Detonative  Combustion  of  Spherical  and  Non-spherical  Source  Regions 

The  work  reported  in  this  section  has  resulted  in  a  rather  complete 
understanding  of  the  manner  of  which  source  behavior  affects  the  blast 
wave  under  detonative  or  deflagrative  combustion  in  the  source  region. 

The  first  paper  entitled  'The  Blast  Wave  Generated  by  Spherical  Flames," 
by  R.  A.  Strehlow,  R.  T.  Luckritz,  A.  A.  Adamczyk  and  S.  A.  Shimpi , 


6 


Combustion  and  Flame,  ,  pp.  279-310  (1979),  contains  a  complete  description 
of  the  effect  of  a  spherically  expanding  flame  with  various  normal  burning 
velocities  up  to  detonative  combustion  on  the  blast  wave  generated  around 
a  spherical  source  region.  The  effect  of  flamo  acceleration  was  also 
studied  by  arbitrarily  legislating  certain  types  of  acceleration  processes 
and  studying  their  effects.  Additionally,  for  low  velocity  flames,  Taylor's 
theory  for  the  blast  wave  produced  by  a  constant  velocity  spherical  piston 
was  adapted  to  flame  propagation  and  was  found  to  fit  the  numerical 
calculations  quite  well. 

Subsequent  to  that  work  a  principle  first  enunciated  by  Stokes  in 
1849  was  adapted  to  deflagrative  combustion  of  clouds  of  arbitrary  shape. 

This  led  to  a  paper  entitled,  ’’The  Blast  Wave  from  Deflagrative  Explosions; 

An  Acoustic  Approach,"  by  R.  A.  Strehlow,  14th  Loss  Prevention  Symposium, 
AIChE,  p.  145  (1981).  This  paper  shows  that  deflagrative  combustion  of 
clouds  of  large  aspect  ratio  with  edge  ignition  yields  blast  pressures 
which  are  considerably  lower  than  blast  pressures  produced  by  a  flame 
propagating  at  the  same  velocity  as  a  spherical  flame  from  the  center  of 
a  source  region.  Experimental  verification  for  this  effect  was  found  in 
the  open  literature. 

An  extension  of  this  work  on  a  sounder  theoretical  foundation  is 
being  performed  by  M.  S.  Raju  at  the  present  time.  Mr.  Raju  has  developed 
a  rather  simple  hydrocode  which  does  not  give  a  great  deal  of  detail  but 
nevertheless  allows  one  to  calculate  pressures  and  flow  velocities  as 
well  as  shock  wave  positions  for  blasts  from  deflagrative  or  detonative 
combustion  of  clouds  of  arbitrary  shape  (with  cylindrical  geometry).  This 
program  will  lead  to  a  Ph.D  thesis  which  will  be  published  in  the  open 
literature. 

The  above  work  on  deflagrative  and  detonative  explosions  in  the  open 
has  shown  quite  definitely  that  damaging  blast  waves  cannot  be  obtained 
from  ordinary  deflagrative  explosions.  Instead,  one  needs  either  a 
detonative  explosion  or  some  type  of  high  velocity  volumetric  combustion 
process  to  occur  before  a  damaging  blast  wave  will  be  generated  when  a  vapor 
cloud  bui'ns.  This  information,  coupled  with  the  SWACER  mechanism  recently 
discovered  by  John  Lee  and  his  co-woi'kers  at.  McGill ,  leads  one  to  the 
conclusion  that  the  really  damaging  vapor  cloud  explosions  are  due  to 
detonative  behavior  in  the  cloud. 
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In  addition  to  the  above  work,  a  paper  entitled,  "Accidental  Explosions," 
was  published  by  Roger  A.  Strehlow  in  the  American  Scientist  magazine, 
July-August,  (1980).  This  paper  represents  a  generalized  summary  of  the 
behavior  of  all  types  of  accidental  explosions  and  a  logical  characterization 
of  them  into  9  catagories. 


Modeling  of  DDT  in  Porous  High  Energy  Propellants 


The  reactive  two-phase  flow  modeling  to  study  the  transition 
potential  of  high  energy,  but  porous  solid  propellant  to  a  steady  detonation 
wave  has  for  the  first  time  shown  that  a  detonation  solution  is  possible. 

This  represents  a  major  achievement,  since  until  this  year,  no  such 
solutions  were  ever  shown.  Improvements  in  the  analysis  which  contributed 
to  the  successful  DDT  modeling  included: 

(a)  several  significant  advances  in  the  numerical  integration  techniques 
to  handle  extreme  shock  wave  formation; 

(b)  use  of  more  appropriate  dynamic  gas  permeability  constitutive 
relations,  based  on  AFOSR  supported  research  under  this  grant; 

(c)  utilization  of  improved  shock  Hugioniot  data  for  HMX  detonation 
products,  thereby  modifying  both  the  non-ideal  gas  phase  and 
solid  phase  equations  of  state. 


Appendix  A  is  a  copy  of  a  paper  (to  appear  in  Combustion  and  Flame) 
which  summarizes  our  results  fur  the  DDT  modeling  effort.  The  paper  was 
condensed  from  a  recent  AFOSR  Interim  Report  (Univ.  of  Illinois  Tech.  Report 
AAE  81-1),  by  P.  Barry  Butler  and  H.  Krier,  entitled,  "Shock  Development 
and  Transition  to  Detonation  Initiated  by  Burning  in  Porous  Propellant 
Beds."  The  paper  presents  solutions  which  indicate  clear  regimes  for 
detonation  potential.  Predictions  for  the  run-up  length  to  detonation  are 
shown  bounded  as  a  function  of  either  propellant  chemical  energy  content, 
propellant  burning  rate,  bed  porosity,  or  particle  size  (granulation). 
Comparison  of  DDT  with  limited  data  available  in  the  literature  indicates 
gcod  qualitative  agreement  with  many  of  these  predictions  and  exceptional 1> 
good  quantitative  agreement  with  detonation  pressure  and  detonation  speeds, 
as  a  function  of  propellant  bulk  density  and  propellant  energy. 
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F .  Measurements  of  Gas  Permeability  in  Packed  Beds 

This  past  year  additional  experiments  have  been  carried  out  that 
extend  the  range-of-validity  of  the  packed-bed  gas-permeability  friction 
coefficient,  a  relation  of  prime  use  in  the  gas  dynamics  associated  with 
the  DDT  modeling.  This  correlation  was  tested  on  both  data  obtained  in  our 
laboratory,  as  well  as  data  (downstream  pressure  decrease  in  packed  beds 
of  solid  particles)  from  other  laboratories. 

The  gas  permeability  correlation  matches  the  relatively  low  Reynolds 
number  conditions  of  published  work  in  the  literature,  but  clearly 
indicates  that  one  should  not  extend  those  conventionally  used  correlations 
to  relatively  high  Reynolds  number  conditions.  Basically,  our  correla¬ 
tion  indicates  that  the  bed  is  significantly  more  permeable  than  would  be 
predicted  by  using  classical  correlations  (generally  found  in  the  chemical 
engineering  literature).  A  technical  paper  is  currertly  being  written 
for  submission  to  a  journal  in  the  open  literature  which  will  summarize 
our  findings  in  this  area. 

Dynamic  measurements  of  high  pressure  decay  through  packed  beds  were 
also  made  using  a  blow-down  air  facility,  so  that  the  initial  pressure 
decayed  from  2500  psia  (17.2  MPa)  to  100  psia  (0.69  MPa)  in  a  matter 
of  several  seconds.  The  transient  test  data  indicate  that  our  steady- 
state  gas  permeability  correlation  was  still  valid  for  this  type  of  dynamic 
event.  Although  transient  gas  permeability  data  in  sub-millisecond  times 
would  best  confirm  to  our  DDT  flow  modeling  regime,  we  believe  that  until 
such  tests  are  made  the  gas-particle  drag  correlation  now  in  use  represents 
the  best  state-of-the-art  relation. 
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NEW  DISCOVERIES 


A.  Reactive  Gas-Air  Mixtures 

1.  Direct  initiation  of  detonation  in  numerical  simulation  can  be 
achieved  by  proper  distribution  of  a  chemical  accelerator  in  a  reactive 
fuel-air  cloud. 

?.  Limiting  values  of  the  controlling  parameters  in  the  chemical 
accelerator  model  exist. 

3.  Source  behavior  affects  the  possibility  of  transition  to  detonation. 

4.  Either  large  flame  area  rate  increases  or  flame  accelerations  are 
needed  to  cause  transition  to  detonation. 

5.  Edge  ignition  provides  so  much  lateral  relief  that  relatively 
low  blast  pressures  are  produced. 

6.  Ordinary  deflagrative  sources  should  not  produce  severely- 
damaging  blast  waves. 

B.  High  Energy  Propellants 

1.  For  the  first  time  complete  deflagration  wave  to  detonation  wave 
transition  has  been  shown  as  a  solution  to  the  reactive  gas-solid 
dynamics  model  of  the  flow  physics. 

2.  Detonation  solutions  require  critical  values  of  chemical  energy, 
bed  porosity,  propellant  pyrolysis  rates,  and  particle  granulation. 

3.  Solid  particle  compressibility  is  a  requirement  in  order  to 
achieve  DDT  in  moderate  porosity  propellant. 

4.  The  burning-rate  pressure  index  of  high  energy  propellants  should 
not  exceed  values  of  3/4  for  HMX  energy  levels  and  average  particle  sizes 
greater  than  250  pm  in  diameter.  Detonation  is  not  predicted  for  indices, 
n,  less  than  0.80,  for  cases  studied  to  date. 

5.  For  high  energy  propellants,  packed  to  porosities,  $  ,  ranging 
from  0.25  to  0.50,  no  detonation  transition  is  predicted  if  the  average 
particle  size  is  greater  than  500  pm  (4  mm).  Hence,  friable  propellants 
should  be  avoided  to  preclude  DDT  hazard. 
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APPENDIX  A: 

MODELING  OF  SHOCK  DEVELOPMENT 
AND  TRANSITION  TO  DETONATION  INITIATED 
BY  BURNING  IN  POROUS  PROPELLANT  BEDS 

by 

H.  Krier,  P.  B.  Butler,  M.  F.  Lembeck 
Department  of  Mechanical  and  Industrial  Engineering 
University  of  Illinois  at  Urbana-Champaign 


ABSTRACT 

This  paper  deals  with  the  analyses  of  deflagration  to  detonation 
transition  (DDT)  occurring  in  a  packed  bed  of  granular,  high  energy  solid 
propellant.  A  reactive  two-phase  flow  model  of  this  phenomena  is  solved 
by  utilizing  a  Lax-Wendroff  finite  differencing  technique.  Utilizing  an 
appropriate  gas  phase  nonideal  equation  of  state  and  high  pressure  gas 
permeability  relations  with  an  improved  numerical  integration  technique, 
one  can  predict  the  transition  to  a  steady  detonation  from  initiation  by 
deflagration. 

Analyses  are  presented  that  clearly  indicate  the  effect  of  the 
propellant  physical  and  chemical  parameters  on  the  predicted  run-up  length 
to  detonation.  Predictions  of  this  run-up  length  to  detonation  are  presented 
as  a  function  of  propellant  chemical  energy,  burning  rate,  bed  porosity, 
and  granulation  (size) .  Limited  comparison  with  actual  DDT  data  in  the 
literature  indicates  qualitative  agreement  with  these  predictions. 


To  appear  in  Combustion  and  Flame  (December  1981) 
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INTRODUCTION 

This  paper  summarizes  the  analysis  associated  with  the  accelerating 
deflagration  wave  in  a  porous  medium  of  reactive  solid  propellant.  The 
hazard  of  DDT  (def lagration-to-detor.ation  transition)  in  solid  propellants, 
especially  solid  propellants  burning  in  rocket  motor  environments,  is 
usually  not  considered  a  possibility.  However,  it  may  be  that  under 
certain  situations  (for  example,  a  grain  structure  failure)  the  solid 
motor  may  crack  and  form  regions  of  granular  or  porous  propellant.  When 
flame  from  the  surface-deflagrating  propellant  reaches  this  seam  of  porous 
material  it  will  accelerate  into  this  medium  and  be  supported  by  convective 
heat  transfer  from  the  burned  gas  into  the  unignited  porous  region.  If,  in 
addition  to  this,  the  product  gases  are  confined  to  a  finite  volume, 
the  accelerating  deflagration  could  transit  into  a  detonation.  Propellants 
exhibiting  a  high  chemical  energy  per  unit  mass  and  capable  of  rapidly 
generating  gases  through  their  burning  rate  are  more  likely  to  experience 
this  type  of  deflagration  to  detonation  transition. 

Analysis  of  the  flows  in  such  an  unsteady  two-phase  mixture  is  a 
complicated  exercise.  Work  has  been  underway  at  the  University  of  Illinois 
since  1976  to  develop  a  reactive  hydrodynamic  code  in  which  the  combus¬ 
tion  of  porous  propellants  can  be  modeled  in  such  a  way  as  to  predict  the 
behavior  of  a  corvectively  driven  flame  in  a  confined  situation.  Details 
of  these  modeling  exercises  are  found  in  Ref.  1-4.  In  an  evolutionary 
manner  this  work  included  the  formulation  of  the  two-phase  flow  conserva¬ 
tion  equations,  first  assuming  that  the  mixture  was  a  continuum,  and  at  a 
later  date  treating  each  phase  as  a  separate  continuum  irrespective  of 
the  mixture  properties.  The  most  recent  analysis  of  the  unsteady  two-phase 
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flows  associated  with, but  prior  to, DDT  is  documented  in  the  paper  by 
Hoffman  and  Krier  [3].  This  work,  therefore,  represents  the  starting 
point  for  the  study  that  is  presented  here. 

The  reader,  after  reviewing  the  above  noted  reference,  can  easily 
see  that  the  modeling  of  this  transient  phenomena  utilizes  a  number  of 
important  constitutive  relations  which  form  closure  of  the  conservation 
equations.  For  example,  one  must  have  information  cn  the  burning  rate 
of  the  material  that  is  a  function  of  the  surrounding  pressure  and  temper¬ 
ature.  One  also  requires  expressions  for  the  dynamic  gas  permeability 
and  the  subsequent  heat  transfer  rates  of  the  hot  gases  as  they  are  forced 
into  the  unignited  porous  material.  In  addition,  relations  which 
represent  the  resistance  to  compaction  of  the  solid  matrix  must  be 
included.  Equations  of  state,  not  only  for  the  high  pressure  in  the  product 
gases,  but  also  for  the  solid  itself,  must  be  supplied  (as  shown  in 
Ref.  3).  The  assumption  of  an  incompressible  solid,  although  providing 
some  reasonable  answers  as  far  as  the  deflagration  speed,  is  not  an 
accurate  indicator  of  the  peak  pressures  that  are  possible  during  the 
accelerating  deflagration  mode.  Since  these  pressures  are  precursors  to 
the  final  detonation  solutions  that  would  be  expected,  it  is  clear  that 
a  compressible  solid  must  be  modeled. 

The  analysis  of  the  DDT  problem  requires  the  solution  of  the  time- 
dependent  conservation  equations  for  both  particle  and  gas  phases  and  the 
necessary  constitutive  relations  noted  above.  The  conservation  equations 
form  a  system  of  nonlinear  hyperbolic  equations,  which  require  numerical 
finite  differencing  schemes.  A  recent  report  has  summarized  several  of 
these  integration  shcemes  and  evaluated  which  are  most  useful  for  this  type 
of  reactive  flow  [5], 
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Previous  Results  on  DDT  Modeling 

A  review  of  Refs.  2-4  indicates  that  a  steady  state  detonation 
solution  was  not  a  predicted  result.  Based  on  the  steady  shock  hydro¬ 
dynamics  for  a  "homogeneous"  reactive  material  with  a  known  chemical 
energy,  initial  solids  loading  and  material  parameters,  a  steady  state 
detonation  (CJ)  would  propagate  at  speeds  of  the  order  of  5-8  mm/ysec, 
with  a  detonation  pressure  of  the  order  of  12-20  GPa  [5].  Although  a 
fairly  rapid  flame  front,  often  approaching  a  steady  state  speed  of 
2  mm/ysec,  was  typical  of  the  solutions  presented  in  the  work  by  Gokhale 
and  Krier[2],  Kezerle  and  Krier  [4],  and  Hoffman  and  Krier  [3],  the 
associated  peak  pressures  were  never  of  significant  manner  to  suggest 
that  these  final  flame  spreading  rates  were  characteristic  of  an  actual 
detonation.  Although  the  peak  pressures  predicted  in  Refs.  2-4  ranged 
between  1/2  to  5  GPa  and  were  consistent  with  such  velocities,  these 
pressures  could  not  represent  the  detonation  (CJ)  pressures. 

Continued  work  (which  is  reported  here)  indicates  that  it  is  now 
possible  to  obtain  a  detonation  solution,  but  in  order  to  do  so  a 
number  of  important  modifications  and  corrections  to  the  work  reported 
in  Refs.  2-4  were  required. 

DDT  in  Two-Phase  Reactive  Flow  (Experimental  Data) 

The  transition  from  deflagration  to  detonation  in  a  porous  reactive 
medium  is  an  unsteady  process.  Hot  gases  generated  from  the  propellant 
surfaces  are  driven  forward  into  theunbumed  solid  matrix  by  the  pressure 
gradient  developed  at  the  ignition  front.  This  phenomena  is  r.ot  found 
when  a  nonporous  solid  detonates,  since  only  pressure  disturbances  can  be 
propagated  ahead  of  the  ignition  front.  It.  is  this  convective  heat  transfer 
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to  the  unignited  propellant  and  the  extended  deflagration  reaction  zone, 
with  the  resultant  compression  of  the  solid, which  makes  DDT  in  a  porous 
bed  unique  when  comparing  it  to  DDT  in  either  an  all  gas  or  an  all  solid 
regime. 

Much  experimental  work  on  DDT  in  porous  material  has  been  carried 
out  in  the  past  decade.  Bernecker  and  Price  have  published  their 
earliest  detailed  results  on  DDT  in  Refs.  6  and  7.  Other  experimental 
studies  prior  to  these  include  the  work  of  Griffiths  and  Grocock  [8]  and 
Taylor  [9] . 

The  work  presented  in  Ref.  7  by  Bernecker  and  Price  is  a  study  of 
DDT  in  RDX  (cyclotrimethylenetrinitramine) ,  a  shock  sensitive  high- 
energy  explosive.  In  their  experiments  the  RDX  was  packed  into  a  thick 
walled  tube  having  an  inside  diameter  of  16  mm  and  being  approximately 
300  mm  (12  in.)  in  length.  Both  ends  of  the  column  were  closed  and 
ionization  probes  were  located  throughout  the  bed  to  trace  the  ignition 
front  locus.  The  RDX  was  packed  in  an  inert  wax  mixture  and  had  a  mean 
particle  size  of  200  urn  in  diameter. 

Their  experimental  work  showed  a  convective  ignition  front  traveling 
at  subsonic  velocities  (0.3-0. 9  mm/ysec''  for  most  cases  when  the  initial 
porosities,  4>q,  ranged  between  values  of  0.1  to  0.3  and  in  cases  where 
transition  occurred,  detonation  velocities  ranged  from  6-S  mm/ysec.  The 
length  of  porous  propellant  required  for  compressional  waves  to  overtake 
the  ignition  front  and  transit  into  a  detonation  was  defined  as  2^,  the 
run-up  length.  Typical  values  for  2^  ranged  from  100  to  200  mm. 
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THIi  MODEL  AND  ANALYSTS 


Introduction 

The  analysis  that  follows  attempts  to  model  the  situation  in  which 
a  bed  of  tightly  packed  granular  propellant  is  ignited  at  one  end.  Both 
ends  of  the  packed  bed  are  considered  closed,  thus  treating  the  problem 
when  the  velocity  of  the  rear  wall  is  much  less  than  the  CJ  product 
velocity.  For  this  case,  after  DDT,  the  detonation  propagates  through 
the  unignited  region  at  the  CJ  detonation  velocity  and  detonation  pressure 
and  is  followed  by  an  isentropic  expansion  of  the  gases  to  the  pressure 
at  the  stationary  wall.  The  gas  surrounding  the  particles  at  the  initial 
time  is  considered  inert,  and  to  be  at  atmospheric  pressure.  It  is  also 
assumed  that  the  inert  gas  will  fully  mix  with  the  gases  being  generated 
from  combustion  of  the  pr  pellant  as  time 'progresses. 

For  numerical  simplicity  the  propellant  particles  are  assumed  to  be 
unisized  spheres.  Particles  of  interest  range  in  diameter  SOum  < 'd  <  200um. 
To  treat  multi-sized  particles,  one  would  require  N  independent  equations 
of  mass,  momentum  anti  energy  for  the  solid  phase,  where  N_  is  the  number 
of  initially  different-sized  particles.  A  solids  loading  of  74%  is  the 
tightest  possible  for  unisized  spheres,  obtainable  by  arranging  the  spheres 
in  a  face  centered  cube.  However,  assuming  granular  deformation  occurs 
under  high  stress  loads,  greater  solids  loadings  may  be  predicted  without 
error.  Obviously,  the  spherical  geometry  must  be  altered  for  this  to  occur. 

As  the  small  fraction  of  propellant  particles  ignited  at  time  t=0 
burn,  hot  gases  are  generated  as  a  function  of  the  pressure-dependent 
burning  rate  law  and  surface-to -volume  ratio  (3/r  )  of  the  spheres. 

These  hot  gases  are  convected  forward  through  the  lattice  of  unburned 
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propellant  and  flow  gradients  develop,  as  dictated  by  the  solution  of  the 
conservation  equations  and  the  necessary  constitutive  relations. 

Heat  transfer  from  the  hot  gases  to  the  unignited  propellant  particles, 
dependent  on  the  velocity  of  the  gas  relative  to  the  particles  and  several 
gas  properties  (i.e.  viscosity,  thermal  conductivity),  transports  energy 
from  the  gas  to  the  solid  phase.  Subsequent  ignition  of  particles 
further  down  the  bed  is  assumed  to  occur  when  a  critical  solid  phase 
internal  energy  is  reached  {2].  This  energy  can  be  expressed  as  a  critical 
increase  in  solid  phase  temperature ,  T  ,  since  the  specific  heat  of  the 
solid  is  assumed  to  be  kno<m. 

As  time  progresses  the  gas  pressure  behind  the  ignition  front  increases 
due  to  the  confinement  of  the  gases  from  the  closed  rear  boundary  and  the 
pressure-dependent  rate  of  mass  generation  in  the  gas  phase.  Under  certain 
conditions*  the  pressure  gradient  can  develop  into  a  shock  front  which 
overtakes  the  ignition  front  propagating  through  the  bed.  When  this  occurs 
the  ignition  front  experiences  the  transition  from  deflagration  to  detonation. 

At  the  transition  point  the  ignited  region  (zone  of  gas  generation) 
narrows  in  width  and  is  followed  by  a  region  of  all  gas  where  the  propellant 
particles  are  completely  burned  out,  as  depicted  in  Figure  1.  The  thick¬ 
ness  of  the  reaction  zone  is  a  function  of  the  initial  particle  size  and 
solids  loading,  and  thicknesses  approaching  1mm  may  be  possible  as 

r  -*■  0  and  d>  •>  0. 
p  ro 


*  These  conditions  will  depend  upon  the  solid  chemical  energy,  granulation 
(size  and  loading),  burning  rate,  ignition  energy,  etc. 


Assumpt ions 

In  order  to  numerically  model  DDT  in  two -phase  flow,  while  retaining 
the  physics  of  the  problem,  several  key  assumptions  had  to  be  made. 

There  assumptions  are  similar  to  those  made  by  previous  investigators 
(Refs  2  and  4): 

(1)  Both  the  solid  and  gas  phases  are  indepen<  ently  treated  as  con¬ 
tinuums  requiring  their  own  conservation  relations. 

(2)  Each  phase  interacts  with  the  other.  This  is  modeled  by  the 
mass,  momentum  and  energy  interaction  terms  in  the  conservation 
equations. 

(3)  All  propellant  particles  are  unisized  spheres. 

(4)  Ignition  of  a  propellant  particle  is  obtained  when  a  critical 
energy,  expressed  as  a  particle  temperature,  is  transferred  to 
the  solid. 

(5)  The  propellant  particles  are  initially  surrounded  by  an  inert 

gas  at  temperature,  T  . 

®o 

(6)  During  combustion  of  the  propellant,  the  gaseous  products  mix 
with  the  inert  gas  described  in  assumption  #5. 

(7)  Both  ends  of  the  bed  are  closed  allowing  no  gases  to  escape. 

(8)  The  specific  heats  at  constant  volume,  Cv>  for  both  phases  are 
constant . 

(9)  When  the  solid  phase,  at  a  given  x-location  in  the  bed,  displays 
a  porosity  <p  >  0.98  and  a  particle  volume  less  than  one-tenth  of 
its  initial  volume,  it  is  burned  out  and  no  longer  generating 
gas.  Results  show  this  phenomena  to  proceed  smoothly  from  the 
left  boundary  and  thus  not  leaving  any  'holes'  in  the  continuum. 
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This  assumption  was  necessary  to  prevent  a  singularity  from 

arising  as  r  -*•0  and  T  ■*  *. 

P 

(10)  All  the  product  gases  obey  an  assumed  nonideal  equation  of 
state,  which  is  reviewed  in  detail  below. 

(11)  The  solid  particles  are  compressible,  without  heating  up,  obey¬ 
ing  a  modified  Tait  equation  of  state. 

(12)  Once  ignited,  the  particles  are  assumed  to  burn  on  the  outer 
surface  only,  at  a  known  pressure-dependent  rate  law. 

(13)  At  some  initial  time,  a  "narrow"  region  at  one  end  is  ignited, 
burning  at  the  low  pressure  prescribed. 


Governing  Equations 

Numerical  modeling  of  the  two-phase  reactive  flow  process  in  DDT 
involves  the  conservation  of  mass,  momentum  and  energy  per  unit  volume 
in  both  gas  and  solid  phases.  This  is  a  system  of  six  conservation  condi¬ 
tions  which  form  a  set  of  nonlinear  hyperbolic  partial  differential 
equations  coupled  by  the  interphase  mass,  momentum  and  heat  transfer 
terms  (r,  P,  Q).  The  conservation  equations  for  two-phase  reactive  flow 
have  been  developed  previously  and  Refs.  3  and  4  will  provide  the  reader 
with  the  definitions,  assumptions  and  expressions  for  the  six  following 
field-balanced  conservation  equations.  Also,  Appendix  A  will  provide  an 
explanation  of  the  interphase  transfer  terms. 


Gas  Continuity 
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Particle  Continuity 
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The  subscripts  g  and  p  denote  gas  and  particle  respectively.  In  Eqs.  1  to 
6,  the  phase  densities,  and  ,  are  defined  as 

P1  =  pg4>  and  P2  =  (l-<j>-)p  (8) 

The  porosity,  <p,  is  defined  as  the  ratio  of  the  instantaneous  gas  volume 

to  the  mixture  volume.  Hence,  the  solids  fraction  is  (l-<j>). 

In  addition  to  the  six  conservation  equations,  three  constitutive 

relations  are  needed  in  order  to  solve  for  the  nine  unknown  variables; 

p  ,  p  ,  u  ,  u  ,  T  ,  T  ,  <p,  P  and  P  ,  These  relations  include  state 
g  P  g  P  g.P  •  '  g  P 

equations  for  both  gas  and  solid  states  and  a  stress-resistance  relation 
for  Pp.  Appendix  A  gives  a  complete  listing  of  the  relations  use. 


Improvements  in  Modeling 

Since  the  work  reported  in  Ref.  3,  certain  "improvements"  in  the 
modeling  effort  have  allowed  solutions  which  may  be  considered  to  be 
actual  detonations.  These  improvements  are  discussed  in  some  detail 
later  in  the  text,  but  basically  include: 

1.  Implementation  of  the  necessary  gas  phase  (nonideal)  equation 
of  state,  to  insure  that  at  the  CJ  (Chapman-Jouget)  conditions, 
the  isentrope  provides  for  a  "gamma  law"  suitable  at  the  hydro- 
dynamic  CJ  state. 

2.  Implementation  of  a  new  gas-particle  friction  coefficient,  as 
developed  by  Wilcox  and  Krier  [10],  for  flows  at  the  high  Reynolds 
numbers  encountered  in  the  developing  DDT  flows.  Previously  such 
coefficients  were  based  on  data  only  available  for  moderate  Reynolds 
number  ranges. 
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3.  Utilization  of  an  improved  numerical  integration  scheme,  which 
allowed  for  a  reduction  in  the  grid  spacing  (and  hence  reduction 
of  the  time  increment)  without  the  usual  penalty  of  excessive 
computation  costs  and  numerical  instability  that  often  follows 
when  the  total  number  of  integrations  is  significantly  increased. 
Details  are  documented  in  a  recent  report  by  Butler  and  Krier 

[Sj. 

Equations  of  State:  Constraint  Due  to  the  Detonation  State 

As  has  been  mentioned,  a  nonideal  equation  of  state  must  be  utilized 
for  the  product  gases.  The  analysis  presented  here  uses  a  nonideal 
equation  of  state  for  hard  spheres  suggested  by  S.  J.  Jacobs  [11].  Pre¬ 
vious  to  this  study,  a  covolume-type  state  equation  with  data  made  available 
by  Cook  [12]  was  used  (see  discussion  in  Ref.  2). 

The  hard  sphere  equation  of  state  takes  the  form 

j~r  =  1.0  +  bp  +  c  (bp)2  +  ...  (9a) 

where  the  constants  b  and  c  are  determined  by  the  value  of  the  gamma  law 
coefficient,  y,  for  the  product  gases.  Here,  y  is  the  negative  logarithmic 
slope  of  the  isentrope  tangent  to  both  the  Rayleigh  and  Hugoniot  lines 
at  the  CJ  point  in  the  detonation  state.  That  is,  the  slope  of  the  isentrope 
that  the  gases  expand  along  in  the  product  state.  The  reader  may  refer 
to  Ref.  5  for  the  complete  solution  of  the  constants  b  and  c. 

Values  for  y  range  from  two  to  three  for  detonating  high  density 
explosives  [13].  For  the  baseline  case  considered  in  this  study,  a 
value  of  y  =  2.05  was  selected  and  the  corresponding  nonideal  equation  of 
state  is: 


i 
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When  the  above  coefficients  (2.S  and  0,50)  were  altered  to  treat  a  case 
for  y  =  3.0,  excessively  high  gas  and  particle  temperatures  were  predicted, 
as  one  might  expect.  In  addition,  during  the  numerical  integration,  severe 
oscillations  in  gas  and  particle  temperatures  occurred  in  most  cases 
when  Y>  3.  Since  one  is  always  constrained  by  the  numerical  integration 
schemes  that  are  employed  to  handle  very  severe  gradients  in  the  flow, 
it  is  understandable  that  previous  efforts  in  DDT  modeling  [2,  3],  which 
were  utilizing  such  high  y  values  for  the  product  gases,  almost,  always 
ran  into  numerical  integration  problems.  For  example,  evaluating  the 
covolume  state  equation  previously  used  in  Ref.  3,  a  value  of  y  =  3.6 
is  calculated.  This  large  value  for  y  may  be  one  important  reason  why 
the  calculations,  as  reported  in  Refs.  2-4,  were  unable  to  handle  the 
high  pressures  associated  with  steady  state  detonations. 

In  the  solid  phase  the  particles  obey  a  modified  Tait  equation  allow¬ 
ing  for  compression  of  the  granules.  This  is  written  as: 


P_  =  ^ 
P  P, 


3P 

P  + 
K 

•  o 


1/3 


(10) 


where  Kq  is  the  bulk  modulus.  Ref.  3  discusses  the  Tait  equation  in  more 

5 

detail.  A  typical  value  is  Kq  =  1.38  GPa  (2.0  x  10  psi). 


CJ  (Detonation)  State 

In  the  text  by  Fickett  and  Davis  [13]  equations  are  presented  for 
estimating  detonation  velocity,  D^j ,  and  the  detonation  pressure,  P^j,  for 
a  steady  state  detonation.  These  are  given  for  a  single  phase  explosive: 
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Here,  p  is  the  initial  solid  material  density  and  E„,ir.,,  is  the  chemical 
P  CHEM 

o 

energy  liberated  by  burning  the  solid.  Since  the  problem,  being  considered 

is  two-phase  (solid-gas),  Eqs.  11  and  12  must  be  modified  to  account  for 

this  by  converting  p  ,  the  initial  solid  density,  to  p-  ,  the  initial 

*o  /  o 

solid  phase  density  p^  =  p  (l-ijj^).  For  an  initial  solid  density  p  =  1994 

o  Po  0  Po 

3 

kg/m  ,  an  initial  porosity  <f)Q  =  0.30,  a  chemical  energy  =  5.48  MJ/kg, 

and  assuming  y=  2,05,  Eqs.  11  and  12  give  respectively:  =  16.06  GPa 

and  D_T  =  5.92  . 

CJ  ysec 

These  equations  were  developed  from  the  jump  equations  for  one-phase 
flow  where  the  equation  of  state  for  the  product  gases  was  assumed  ideal 
Because  of  this,  the  above  values  should  be  considered  only  as  good  esti¬ 
mates  for  the  detonation  pressure  and  velocity. 


Gas  Permeability 

One  of  the  key  constitutive  relations  required  in  the  analysis  is 
the  gas-particle  (interphase)  viscous  force,  which  governs  hot  gas  pene¬ 
tration  into  the  unignited  region  of  the  granular  material.  As  presented 


in  Appendix  A: 


(u  -  u  )  f 
2  v  g  p  pg 


where  f  is  the  drag  coefricient .  Until  recently,  the  packed  bed  corre- 

pg 

lations  by  Ergun  or  Kuo  and  N’ydegger  (as  reviewed  in  Ref.  10)  were  uti¬ 
lized  for  f  .  Thus,  the  modeling  efforts  presented  in  Ref.  5  and  4  used 
Pg 
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the  expression  of  Kuo  and  Nydegger  [15]: 

2  87 

fpg  ■  5  )  04) 

Eq.  14  was  developed  for  460  <  Re  <  14,600.  Here,  Re  is  the  appropriate 
Reynolds  number,  defined  as: 

Re  -  [(4>ug)Pg  2rp] /Pg  (15) 

Based  on  experiments  at  both  high  gas  velocities  and  high  Reynolds  .'■•umbers. 
Wilcox  and  Krier  [10]  developed  the  correlation: 


f  =  5.06  x  105  r  Re/u2 
Pg  P  g 


(16) 
,5 


where  u  is  given  in  m/sec,  r  in  meters,  and  the  constant  5.06  x  10" 
g  P 

2 

must  have  units  of  m/s  .  This  expression  is  not  valid  for  either  very 


low  gas  velocities,  since  Eq.  16  would  give  f  -*■  00  as  u  -*•  0,  or  for 

r&  o 


very  high  gas  velocities*,  since  f  0  as  u  +  ®.  The  equation  was 

pg  g 

found  to  be  fairly  accurate  for  10^  <  Re  <  (2  x  10*’)  and  15  m/s<u  <150  m/s 

o 


While  a  straight  forward  comparison  is  not  easy,  the  difference  between 


the  value  for  f  as  predicted  by  Eq.  14  versus  that  by  Eq.  16  can  be  seen 
in  the  example.  At  Re  =  10^  and  $  =  0.4,  Eq.  14  gives  f  =  55,000  and 

r  o 


Eq.  16  gives  a  value  of  f  =  5448.  To  utilize  Eq.  16  one  must  specify 


the  average  gas  velocity,  u  ,  and  the  particle  radius,  r  ,  that,  were  used 

g  P 

to  obtain  the  Reynolds  number.  A  kinematic  viscosity  u  / p  =  1.8  x  10 

6  o 


m  /sec,  a  particle  radius  r  =1.0  mm,  and  an  average  gas  velocity 

4 

u  =  30.5  m/sec  were  used  to  obtain  a  Reynolds  number  Re  =  10  .  From  this 
g 


*Data  from  Ref.  10  was  limited  to  u  <  300  m/sec. 
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particular  example,  the  Wilcox/Krier  correlation  (Eq.  16)  allows  about  10 
times  the  permeability,  i.e.,  1/10  the  viscous  drag  force  as  correlated 
by  the  Kuo/Nydegger  relation. 

Solutions  to  the  flow  process  leading  to  DDT,  discussed  in  the  follow¬ 
ing  chapter,  clearly  indicate  that  sufficient  gas  permeability  is  necessary 
to  allow  for  a  detonation  transition. 


NUMERICAL  INTEGRATION 


Finite  Difference  Mesh 

To  solve  the  Eulerian  formulated  system  of  conservation  equations 

discussed  above  with  the  constitutive  relations  listed  in  Appendix  A,  the 

length  of  the  bed  being  integrated  over  is  divided  into  I  segments, 

each  a  const  ait  Ax  in  width  (i.e.,  x^  =  jAx;  j  =  1,  2,  3. ...I). 

At  time  t  =  0  values  of  the  nine  independent  variables;  p  ,  p  , 

2  p 

u  ,  u  ,  T  ,  T  .  P  ,  P  and  d>  are  initialized  at  each  jth  x-location  in 
g  P  g  P  g  P 

the  grid.  Before  incrementing  the  primary  variables  (i.e.  mass,  momentum 
and  energy)  to  the  future  time,  t  =  t  +  At,  the  auxiliary  variables  in 
the  equations  (e.g.,  drag,  gas  generation,  heat  transfer)  must  be  computed 
at  the  present  time,  t  =  tQ.  The  nine  equations  are  then  solved  at  the 
incremental  time,  t  =  t  +  At,  by  a  modified  Lax-Wendroff  finite  differenc¬ 
ing  scheme.  This  method, along  with  another  used,  are  presented  in  Ref. 

5.  The  second  method,  developed  by  Rubin  and  Berstein  [16],  was  implemented 
for  several  test  cases  and  gave  results  similar  to  that  of  the  Lax- 
Wendroff  scheme. 

The  time  increment.  At,  over  which  the  equations  are  solved  is 
calculated  by  the  Courant ,  Fredrichs,  Levy  stability  criteria,  for 


hyperbolic  equations,  i.e.. 


=  ■fe+FI)  (17) 

In  Eq.  17  the  term  c_  is  the  mixture  sound  speed  and  |u|  is  the  maximum  gas 
velocity  in  the  bed.  Also,  A  is  a  stability  constant,  less  than  unity. 

For  most  cases  A  =  0.5  was  used.  This  smaller  time  increment  was  necessary 
in  order  to  integrate  the  equations  when  large  gradients  developed  in 
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Initial,  and  Boundary  Conditions 

To  initialize  the  problem  the  bed  is  assumed  to  be  quiescent,  i.e., 
at  a  constant  gas  temperature  and  constant  gas  pressure  throughout  the 
length  of  the  bed.  The  spherical  propellant  particles  are  typically  fixed 
at  a  constant  solids  loadings,  although  a  variable  initial  porosity  can 
be  treated.  Then  to  initiate  the  flow,  the  propellant  at  the  first  few 
grid  points  is  assumed  ignited  and  generating  gas.  To  be  consistent 
with  the  flat  initial  pressure  profile,  all  gas  and  particle  velocities 
at  time  t  =  0  are  set  equal  to  zero.  Table  1  is  a  summary  of  typical 
input  data  for  the  cases  studied. 

Modification  to  the  Integration  Scheme 

To  solve  the  flow  on  a  digital  computer  (which  is  to  represent  a  DDT 
phenomena)  requires  the  repetitive  integration  of  many  equations.  All 
of  the  conservation  equations  and  constitutive  relations  are  solved 
individually  at  each  grid  point  in  the  entire  bed,  for  each  time  increment. 

A  study  of  several  test  cases  showed,  for  instance,  that  when  the 
j-r.ition  front  was  at  a  given  x-location,  there  was  little  activity  several 
grid  points  ahead  of  it.  That  is,  the  gas  and  particles  were  at  velocities 
close  to  zero  and  all  transport  coefficients  were  negligable  in  the  region 
not  far  in  front  of  the  ignition  front.  This  phenomena  occurred  for  both 
the  detonation  state  where  the  front  was  propagating  between  5-8mm/psec, 
and  the  deflagration  state  where  the  propagation  velocity  was  much  less. 
Because  of  this,  an  addition  was  made  to  the  computer  code  which  allows 
the  code  to  only  integrate  the  active  region  and  bypass  the  inactive 


zone . 
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Ahead  of  the  ignition  zone,  the  computer  code  located  the  nearest 
point  to  the  zone  where  there  was  no  significant  particle  or  gas  movement. 

This  point  was  then  designated  p s  the  new  front  boundary  of  the  integration 
region  for  that  particular  integration  step.  When  the  pressure  gradient 
increased,  the  ignition  front  moved  rather  rapidly  through  the  bed  (5-10  ~ — ■) 

(JwCv 

and  a  new  integration  boundary  had  to  be  located  after  each  time  increment. 

Comparisons  of  this  modified  integration  technique  with  a  standard 
fixed  domain  integration  are  reported  in  Ref.  5  and  clearly  indicate  that 
the  predicted  results  are  almost  identical.  To  assure  a  correct  solution, 
the  actual  front  integration  boundary  was  extended  a  few  grid  points  beyond 
the  location  calculated  by  the  code.  The  addition  of  this  new  logic  to 
the  current  code  reduced  computation  time  by  at  least  one  half,  for  the 
same  time  increment.  At,  and  grid  space,  Ax. 

Artificial  Smoothing 

Inherent  to  the  solution  of  the  system  of  interdependent  conservation 
equations  is  a  numerical  instability.  A  small  perturbation  can  in  some 
cases  amplify  with  each  time  increment  and  eventually  destroy  the  numeri¬ 
cal  solution.  This  phenomena  can  start  at  the  first  time  increment  and 
in  ten  to  twenty  integrations  the  oscillations  can  be  so  large  that  the 
solution  becomes  unstable  and  terminates.  In  order  to  smooth  out  these 
oscillations  before  they  amplify,  an  artificial  smoothing  routine  was 
incorporated  into  the  code. 

From  experience  in  integrating  the  two-phase  flow  equations  (Eqs.  1-6) 
with  significant  nonlinear  source-sink  terms,  the  problem  of  numerical 
instability  occurs  often  enough  to  warrant  artifical  smoothing. 

The  analyzation  of  numerous  test  cases  has  shown,  suprisingly,  that 
the  stable  solution  did  not  require  smoothing  of  all  the  variables.  It 
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is  obvious  that  the  following  results  have  an  inherent  dependency  on 
smoothing.  Extreme  care  has  been  taken  to  minimize  these  effects  on  the 
qualitative  trends  and  quantitative  results  predicted.  Nevertheless,  it 
should  be  obvious  that  smoothing  techniques  can  supply  variability  in  the 
predicted  parameters  which  are  reflections  of  the  scheme  utilized,  and 
not  necessarily  of  the  conservation  equations . 


I 
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RESULTS  COMPUTED 


Introduction 

This  section  will  present  the  calculations  made  for  possible  DDT 
occurring  in  packed  beds  of  high  energy,  granulated,  unisized  propellants 
or  explosives.  It  is  obvious  that  there  are  a  very  large  number  of 
loading  combinations  possible  for  our  DDT  study.  Fortunately,  the  work 
of  Hoffman  and  Krier  [3]  and  Krier  and  Gokhale  [2],  in  which  conditions 
of  rapid  convective  flame  spreading  have  been  calculated,  is  available, 
and  can  be  used  as  a  starting  point.  (It  should  be  noted  that  none  of 
the  calculations  made  in  Ref.  5,  2  or  3  predicted  POT).  As  pointed  out  befor 
there  are  probable  reasons  why  this  has  not  been  accomplished  and  it  is 
expected  that  those  improvements  and  modifications  discussed  in  the  previous 
section  will  now  allow  for  the  calculation  of  the  steady  state  detonation 
solution. 

The  study  made  of  a  potential  DDT  attempts  to  model  a  long  column  of 
granulated  material  in  a  closed  pipe  ignited  by  an  energetic  ignition 
material  at  one  end.  In  order  to  model  the  flow  transient,  one  must 
assure  that  the  length  of  the  bed  exceeds  the  run-up  length  to 

detonation.  Since  the  experimental  work  of  Bernecker  and  Price  [6] 
indicated  that  a  10  in  (2S  cm)  bed  was  sufficient  for  most  of  their 
experiments  where  DDT  occurred,  this  length  was  included  as  a  condition 
to  be  treated. 

Since  unisized  spheres  are  bning  treated,  the  initial  porosity  can 
be  no  less  than  =  0.26,  although  randomly  packed  unisized  spheres 
generally  give  a  high  porosity,  about.  t>Q  =  0.40.  Therefore,  for  this 
study  0.26  <  ^  <  0.40.  The  initial  particle  radii  studied  were  also  in 
the  same  range  as  those  considered  by  Refs.  2  and  3.  Results  in  these 


studies  showed  that  particles  must  be  less  than  one  millimeter  in  diameter 
in  order  to  generate  sufficient  gases  for  the  rapid  flame  spreading 
phenomena. 

The  chemical  energy  of  the  material  considered  is  in  the  range  of 
explosives  or  high  energy  propellants  of  the  nitramine  family  (i.e.,  HMX, 
RDX) .  Thus,  the  chemical  energies  studied  were  always  larger  than  4.15 
MJ/kg  (1000  cal/g) .  Other  parameters  one  must  consider  in  the  DDT  studies 
are  the  burning  rate  properties  of  the  propellant.  Again,  the  values  used 
in  Refs.  2  and  3,  which  attempted  to  model  the  burning  rate  of  an  HMX 
solid  propellant,  were  utilized.  However,  the  burning  rate  index,  n, 
is  a  parameter  which  is  explicitly  studied. 

In  this  analysis,  the  deflagration  will  be  initiated  by  assuming  at 
time,  t  *  0,  that  except  for  a  small  portion  of  the  bed  at  one  closed  end, 
the  bed  is  quiescent  and  unreacting.  A  closed  end  situation  is  always 
considered,  and  hence  at  the  two  end  points  (x  *  0,  x  *  L)  it  is  assumed 
that  all  flow  gradients  are  zero  and  that  the  velocities  of  the  particles 
and  gas  must  also  be  zero. 

The  following  section  presents  a  solution  that  indicates  (for  the  first 
time)  a  steady  state  detonation  can  be  predicted.  This  result  is  then 
compared  to  conditions  where  no  such  detonation  solution  occurs.  Additional 


calculations  will  indicate  the  sensitivity  of  the  initial  porosity,  4>q, 
chemical  energy,  burning  rate  index,  n,  and  ignition  energy,  , 

on  the  run-up  length  to  detonation. 


Calculat ions 

Figures  2a  and  2b  present  a  case  where  a  transition  from  deflagration 
to  detonation  has  occurred.  For  this  example,  the  burning  rate  index 
was  n  =  1.0,  the  chemical  energy  =  5.48  MT/kg  and  the  particle  radius 
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r  ■  127  ym  (.005  in) . 

P 

The  pressure-distance  profiles  at  five  separate  times  after  ignition 
of  the  propellant  at  x  ■  0  are  shown  in  Figure  2a.  Examining  the  profile 
for  t  *  SO  usee ,  one  can  observe  that  the  profile  is  characterized  by  a 
shock  front  at  x  «  23  cm  followed  by  a  smooth  expansion  back  to  the 
wall  at  x  ■  0  cm.  The  pressure  in  front  of  the  shock  is  at  atmospheric 
conditions  and,  therefore,  negligible  with  respect  to  the  pressure  behind 
the  shock. 

The  ignition  locus  plot  for  this  particular  case  is  shown  in  Figure 
2b.  Here,  the  ignition  front  moves  through  the  bed  at  a  low  subsonic 
velocity  for  the  initial  ten  microseconds  and  then  accelerates  to  reach 
a  steady  state  velocity  of  DCJ  «  7.2  mm/usec.  This  occurs  within  12  to  15 
cm  from  the  ignited  end.  These  detonation  solutions  for  PCJ  and  DCJ  are 
in  fair  agreement  with  the  approximations  of  Eqs.  11  and  12  for  the  given 

i"Put  <echbc  Po’  V' 

Obviously,  hydrodynamic  steady  state  analysis  (like  Eqs.  11  and  12) 
cannot  guarantee  that  a  transition  from  deflagration  to  detonation  will 
occur.  However,  it  seems  that  "critical"  values  of  porosity  (related  to 
gas  confinement),  gas  generation  rates,  and  chemical  energy  will  provide 
for  a  DDT.  For  example,  when  the  burning  rate  index  was  lowered  to  n  =  0.8, 
as  shown  in  Figures  3a  and  3b,  the  steep  pressure  front  associated  with  a 
detonation  did  not  develop.  Correspondingly,  no  detonation  speed  was 
predicted.  In  this  example  the  peak  pressure  in  the  bed  never  exceeded 
5  GPa,  no  shock  was  predicted,  and  only  a  steady  convection-driven  front 
of  2.2  mm/usec  occurred  at  100  usee  after  ignition  of  x  ■  0. 

According  to  Eq.  11,  for  the  case  where  a  transition  to  detonation 
actually  occurs,  the  steady  state  detonation  pressure,  P^j,  should  increase 
linearly  with  the  chemical  energy,  Also,  Eq.  12  states  that  the 


24 


detonation  velocity,  D^,  is  a  function  of  the  square-root  of  the  chemical 
energy.  Figures  4a  and  4b  present  the  results  of  a  calculation  where  all 
parameters  were  identical  to  those  used  to  give  the  DDT  results  of 
Figures  2a  and  2b,  except  *  6.85  MJ/kg,  an  increase  of  25%.  The 

steady  state  shock  pressure  predicted  for  this  case,  shown  in  Figure  4a 
was  21  GPa.  This  represents  a  (21/16.4  »  1.28)  28%  increase  in  pressure 
over  the  first  case.  The  predicted  detonation  speed  (Fig.  4b),  calculated 
from  the  slope  of  the  x-t  diagram,  was  8.70  mm/ysee.  According  to  Eq.  12, 
the  ratio  of  P  T  for  the  case  given  in  Figure  4b  to  that  presented  in 
Figure  2b  should  be 

/(6.S5/5.48)  »  1,12 

This  is  approximately  the  increase  predicted. 

Table  2  summarizes  the  detonation  pressure,  PCj,  and  the  detonation 
velocity,  D^,  (which  were  the  end  resu!  ■ ....  of  the  DDT  calculations)  all 
as  a  function  of  the  propellant  chemical  energy.  These  predicted  condi¬ 
tions  are  compared  with  the  approximate  steady  state  hydrodynamic  solutions 
discussed  earlier,  i.e.,  Eqs.  11  and  12.  The  excellent  agreement  of  the 
detonation  pressure  with  the  analytic  solution  s.hould  be  noted. 

However,  the  predicted  value  for  the  detonation  velocity,  DCj,  from 
the  hydrodynamic  solution  (Eq.  12)  shows  to  be  slightly  less  than  the  value 
predicted  by  the  code  for  all  cases.  Although  one  cannot  judge  which 
of  the  two  values,  detonation  velocity  or  detonation  pressure,  is  more 
accurate,  the  percent  increase  in  the  detonation  velocity  as  the  chemical 
energy  is  increased  compares  favorably  with  the  hydrodynamic  solution. 
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DDT  Run-up  Length 

The  run-up  length  to  detonation  is  defined  in  this  report  to  be  the 
distance  from  the  closed  end  where  the  bed  is  ignited  to  the  location 
where  both  the  peak  pressure  and  the  detonation  speed  are  constant,  i.e., 
the  equilibrium  steady  state  solution. 

Figures  5  and  6  plot  the  predicted  run-up  length  to  detonation  as 
a  function  of  the  burning  rate  pressure-index,  n  (Fig.  5),  and  initial 
bed  porosity,  <j>Q,  (Fig.  6).  The  "no-solution"  boundary  indicates  that, 
with  the  integration  scheme  used,  the  mesh  si~e  would  have  had  to  be 
drastically  reduced  (thereby  decreasing  the  integration  time  increment) 
in  order  to  obtain  a  stable  solution.  This  is  a  costly  exercise,  but 
future  work  is  planned  to  increase  the  solution  regions.  Figure  5  clearly 
indicates  that,  as  expected,  one  cannot  achieve  detonation  if  the  burning 
rate  during  the  deflagration  phase  is  not  sufficiently  large  (see  "no 
transition"  boundary).  For  the  solids  loading  considered  (3. -^0=  0.70), 
it  would  appear  that  a  minimum  DDT  runup  length  is  5  cm  for  particles 
of  250  ym  in  diameter. 

Figure  6  begins  to  resemble  the  required  "U-shaped"  curve  of 

versus  initial  porosity.  DDT  experiments  by  Korotkov  eb_-aJL-[18)  show 
a  similar  behavior.  One  would  expect  that  for  a  relatively  porous  bed, 

>  0.60,  no  transition  will  occur  since  local  pressure  confinement  is 
limited.  If  the  porosity  is  too  small  no  gas  penetration  for  the  acceler¬ 
ating  deflagration  wave  will  occur.  The  net  result  is  a  porosity  where 
a  minimum  run-up  length  occurs. 

Figure  7  presents  a  study  on  the  effect  of  the  chemical  energy,  > 

on  the  run-up  length,  £_T,  values  which  were  presented  in  Table  2. 

LJ 

As  expected,  the  run-up  length  to  detonation  increases  as  the  amount  of 
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chemical  energy  decreases.  Again,  as  expected  there  is  a  minimum  value 
where  no  transition  occurs,  here  about  =  3.0  MJ/kg.  For  this  case 

it  does  appear  that  a  constant  but  small  run-up  distance  is  still  required 
as  the  chemical  energy  increases  beyond  the  values  studied  here. 


Detonation  Reaction  Zone 

One  measure  of  the  fact  that  detonation  occurs  is  a  plot  of  the 
reaction  zone  width  versus  the  locus  of  the  ignition  front.  This  is  pre¬ 
sented  in  Figure  8  where  the  reaction  zone  is  defined  as  the  region  where 
particles  are  ignited  and  generating  gas  (i.e.,  r  >  0).  As  shown,  the 
zone  initially  increases  during  the  deflagration  phase  as  the  convective 
heat  transfer  provides  energy  to  ignite  more  and  more  of  the  bed.  The 
zone  then  collapses  to  a  thin  (constant)  width  as  the  surrounding  high 
gas  pressure  causes  the  particles  to  burn-out  rapidly.  A  steady  reaction 
zone  thickness  of  approximately  9  mm  is  predicted. 

However,  most  of  the  gas  is  generated  in  a  small  region  immediately 
behind  the  ignition  front  where  the  particles  are  still  relatively  large. 
Note  that  in  all  the  cases  reported  here,  the  initial  particle  diameter 
was  250  ym.  Obviously,  initially  smaller  particles  will  provide  for  a 
thinner  detonation  reaction  zone. 


Comments  and  Interpretations 

The  results  shown  in  Figures  2-8  have  clearly  indicated  that,  as 
expected,  high  solids  loading  of  relatively  small  particle  size  energetic 
propellant  with  confinement  will  transit  into  a  detonation  in  space 
domains  of  several  centimeters.  One  would  expect  that  propellant 
properties  and  packing  configurations  have  limits  where  detonation  cannot 
occur  and  this  was  clearly  shown  in  some  of  the  figures  which  show  the 
run-up  length  versus  property  parameters.  In  conclusion,  it.  will  be  useful 
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to  review  these  studies  to  determine  the  properties  and  configuration.-- 
which  minimize  a  DDT  hazard. 

For  example,  for  a  fixed  chemical  energy,  Eruck.  =  5.48  MJ  kg,  a 

fixed  ignition  energy,  =  KJ/kg,  a  fixed  particle  radius, 

r  =  127  ym  (0.005  in)  and  a  fixed  solids  loading,  (1-0  )  =  0.7,  the 
Po 

burning  rate  index,  n,  must  be  larger  than  n  =  0.86  if  DDT  is  to  occur. 
This  was  shown  in  Figure  5.  Of  course,  had  the  burning  rate  coefficient, 
b,  been  a  different  number,  this  exponent  may  have  been  different.  A 
general  statement  can  then  be  made  that  the  burning  rate,  and  hence  rate 
of  gas  generation,  should  be  kept  as  low  as  possible  to  minimize  DDT. 
Conversely,  the  higher  the  burning  rate,  the  better  the  chance  of  a  transi 
tion  to  detonation  occurring. 

Figure  6  showed  the  run-up  length  to  detonation  versus  initial 
porosity,  <}>o.  As  stated,  for  all  parameters  equal,  there  is  a  maximum 
initial  porosity  where  no  transition  occurs.  Since  a  randomly  packed 
bed  of  unisized  particles  has  an  initial  porosity  of  the  order  <f>Q  s  0.4, 
this  figure  indicates  that  this  is  in  the  region  where  DDT  potential  is 
at  a  maximum.  As  one  reduces  the  initial  porosity  (that  is  increases 
the  solids  loading)  the  run-up  length  to  detonation  slightly  decreases 
until,  as  experimental  work  of  Bernecker  and  Price  [6]  has  indicated, 
there  is  a  minimum  initial  porosity  where  DDT  cannot  occur.  Since  at 
these  initial  solids  loadings  multi- sized  particles  and  mechanical 
packing  are  required,  these  loadings  are  not  of  interest,  whereas  the 
randomly  packed  loadings  are  of  interest. 

A  similar  comparison  of  the  run-up  length  to  detonation  was  shown  in 
Fig.  7  where  the  chemical  energy  content  was  varied.  Recall  that  a  chemi¬ 
cal  energy,  E^^  =  4.18  MJ/kg  (1000  cal/g),  can  be  considered  an  ener¬ 
getic  propellant  material.  Based  on  the  results  shown  in  Fig.  7,  one  can 
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state  that  less  energetic  material  than  this  has  little  chance  of  encounter¬ 
ing  a  DDT.  Doubling  the  chemical  energy  from  F. -  4.18  MJ/kg  to 
^CHEM  "  8.36  MJ/kg,  represents  a  reduction  in  the  run-up  length  of  only 
about  one  half.  To  summarize,  high  energy  propellant  of  the  nitramine 
family, where  >.  4. 18  MJ/kg ,  definitely  fall  within  the  regime  of  a 

DDT  hazard  if  properly  confined. 

The  final  comparison  of  this  type  is  shown  in  Figure  9.  This  shows 

the  run-up  length  to  detonation  versus  the  particle  radius.  As  one  would 

expect,  there  is  a  maximum  particle  radius  (i.e.,  surface-to-volume  ratio) 

where  transition  is  not  predicted  to  occur.  This  is  indicated  by  the 

"no-transition"  boundary  on  that  figure.  The  figure  also  indicates  that 

as  the  particles  get  smaller  in  size,  the  run-up  length  to  detonation 

also  decreases,  as  expected.  Reducing  the  initial  particle  size,  r  , 

^o 

to  values  less  than  25  ym  results  in  gas  generation  rates,  per  unit-volume, 

that  are  so  large  finer  grid  spacing  must  be  utilized  to  assure  stability. 

This  expensive  task  has  been  delayed  and  is  recommended  only  after  improved 

numerical  integration  schemes  have  been  developed. 

A  final  topic  studied  dealt  with  the  effect  of  ignition  temperature 

(or  more  appropriately  AE^^)  on  the  run-up  length  to  detonation.  For 

the  results  shown  in  this  chapter  a  nominal  value  of  T.  =  303°K  was 

ign 

used.  For  this  study  <j>Q  =  0.30,  r  =  127  ym  and  n  =  1.0. 

Calculations  were  made  in  which  T.  varied  over  the  range  294°K  < 

ign  6 

T.  <  350°K.  For  a  constant  initial  bed  temperature  of  T  =  T  =  294°K 
xgn  “  So  Po 

this  represents  a  range  of  ignition  energy  of  0.0  <  AE.  <  25.7  KJ . 


ign  - 


kg 


For  most  of  the  ignition  temperatures  tested  there  was  little  change  in 


the  steady  state  detonation  pressure  or  velocity.  Only  when  AH. 


ign 


approached  the  improbable  value  of  AE^^O  did  the  values  change  significant! > 
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TABLE  1 


Typical  Input  Data 

PARAMETER 
Burning  Rate  Index 

Burning  Rate  Proportionality  Constant 

Initial  Bed  Porosity 
Particle  Diameter 
Bed  Length 
Grid  Spacing 

Chemical  Energy 

Gas  Specific  Heat 

Initial  Bed  Temperature 
Ignition  Temperature  (Bulk) 

Ignition  Energy 


VALUE 


0.8  <  n  <  1.2 

b  *  °-001  fer  ipk>' 


0.25  <  <fr  <  0.50 
—  o  — 

50ym  <_  dp  _<  500ym 
L  =  25.0  cm 


Ax  =  1.27  mm 


4  E  <  7  ^ 
*  kg  ‘them  kg 


rsr  <  c  <  1.9 


kg°k  -  vg  - 


T  =  294°K 
g 

T.  b  303°K 
ign 

E.  b  9.0  ~ 
ign  kg 


"in  ■iiT~r  " 


TABLE  2 


Comparison  of  the  Predicted  Detonation  State  with 
an  Approximate  Hydrodynamic  Solution 


echem 

PT (predicted) 

PCJ(Eq.  11)* 

(predicted) 

DCJ(Eq.  12)** 

0 

CJ 

4.11  MJ/kg 

14.0  GPa 

12.0  GPa 

6.48 

ysec 

5.12  ^ 
ysec 

89  mm 

5.48 

16.6 

16.07 

7.25 

5.92 

51 

6.85 

21.8 

20.08 

8.24 

6.62 

44 

8.22 

24.6 

24.1 

9.17 

7.25 

38 

10.96 

30 . 0 
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10.58 
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2a.  Pressure  history  during  accelerating  deflagration 
in  a  packed  bed  leading  to  a  detonation  transition 
(n  =  1.0) 
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Fig.  2b.  Ignition  front  locus  with  detonation  transition  (n  =  1.0) 
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in  E„irc,,  (compare  to  case  shown  in  Fig.  2a). 

CHEM 


Echem=6.85MJ/KG 
<£o=  0,30 
rp0=  127/xm 
n  =  10 


_Dj  =  8.27 


X  5 


0  10  20  30  40 

Time  (/j. sec) 

Fig.  4b.  Ignition  front  locus  (compare  to  Fig.  2b) 


,  j.  antique 


(cm) 


38 


Porosity  (c pj 


Variation  in  initial  bed  porosity,  i>  ,  on  run-up  length 


No  Transition 


Variation  in  run-up  length  to  detonation,  A  ,  with  chemical 
energy,  ECHBr 


i  »'>.wg|»«mww»,'J  ■  1 . mszm  ,.  1 iwp^_jf)!jiii|p|^^ 


APPENDIX  A 


CONSTITUTIVE  RELATIONS 


As  discussed  earlier,  a  system  of  nine  independent  equations  is 
necessary  in  order  to  solve  for  the  nine  unknown  variables  describing  two- 
phase  flow:  pg>  Pp,  ug>  up,  Eg,  Ep,  Pg,  Pp  and  $. 

Of  the  three  additional  equations  necessary  for  the  solution,  one  is 
the  ninideal  equation  of  state  for  the  gas  phase  as  described  in  the  text. 


p-f-  =  1  +  bpg  +  c  (bpg)2  +  . . .  (A.l) 

g 

Values  for  the  constants  b  and  c  are  discussed  in  Ref.  5. 

The  second  represents  an  equation  of  state  which  relates  the  solid 
density  to  the  stress  on  the  particle.  (See  Ref.  3). 


P 

P 


P  3 


CA.2) 


Finally,  the  third  additional  constraint  is  a  relation  for  the  par¬ 
ticle  phase  stress,  Pp,  as  a  function  of  the  solids  loading  and  the  material 

bulk  modulus,  K  (also  see  Ref.  5). 
o 

As  discussed  in  References  '!  and  3,  one  must  also  specify  functional 
relations  for  the  following: 

A)  T  =  mass  generation  rate  per  unit  volume 


T  =  —  (l-i)  p  r 

r  v  3  p 
P  y 


(A.  o') 


Here,  r  is  the  instantaneous  particle  radius  and  r  Is  t he  surface  burning 

P 

rate  specified  as  a  function  of  pressure  (and  possibly  particle  temperature'. 


rrjJJ.MKius. 


For  all  cases  run  in  this  study: 
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APPENDIX  B 

Accidental  Explosions 


Recent  research  has  significantly  increased  our 
understanding  of  the  causes  and  behavior  of 
accidental  explosions  and  has  produced  new  and 
more  effective  safety  measures 


Accidental  explosions  are  an  unde¬ 
sirable  side  effect  of  man’s  techno¬ 
logical  development.  The  first  such 
explosions  were  undoubtedly  associ¬ 
ated  with  the  manufacture,  handling, 
and  use  of  black  powder.  In  the  Mid¬ 
dle  Ages  alchemists  inadvertently 
caused  explosions  by  mixing  incom¬ 
patible  chemicals.  The  development 
by  Lavoisier  (1789)  of  a  systematic 
nomenclature  for  inorganic  com¬ 
pounds  led  to  the  golden  age  of  pre¬ 
parative  chemistry,  during  which 
many  new  explosive  substances  and 
mixtures  were  prepared,  resulting  in 
more  laboratory-scale  explosions. 
This  in  turn  led  quickly  to  larger  ac¬ 
cidental  explosions  with  the  devel¬ 
opment  of  the  new  field  of  high  ex¬ 
plosive  technology. 

At  about  the  same  time,  the  industrial 
revolution,  with  its  need  for  fossil  fuel 
energy,  introduced  coal  mine  explo¬ 
sions  to  mankind.  The  commerciali¬ 
zation  of  grain  handling  and  milling 
in  the  nineteenth  century  led  to  many 
major  grain  elevator  and  mill  explo¬ 
sions  (Price  and  Brown  1922).  The 
recent  development  of  larger  and 
larger  single-line  petrochemical  pro¬ 
cess  units  and  the  bulk  transport  and 
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storage  of  eve  increasing  quantities 
of  fuels  under ».  <u  pressure  or  at  low 
temperature  have  led  to  boiling-liq¬ 
uid  expanding-vapor  explosions,  or 
BLEVEs,  and  ur.confined  vapor- 
cloud  explosions.  These  can  be  di¬ 
sastrous,  and  their  frequent  incidence 
in  the  last  decade  has  led  to  a  resur¬ 
gence  of  research  into  the  nature  of  all 
explosions. 

The  term  explosion  is  not  very  pre¬ 
cise.  I  shall  use  it  in  this  article  to 
mean  a  process  in  which  the  produc¬ 
tion  or  release  of  a  gas  under  high 
pressure  is  rapid  enough  to  cause  a 
pressure  wave,  usually  called  a  blast 
wave,  to  propagate  through  the  sur¬ 
rounding  medium,  usually  air.  Ex¬ 
plosions  are  audible  at  some  distance 
from  the  source:  nearby  an  explosion 
is  heard  as  a  sharp  “crack,”  while  far 
away,  owing  to  atmospheric  inhomo¬ 
geneities,  it  is  heard  as  a  “boom.” 

The  most  important  aspect  of  any 
accidental  explosion  is  of  course  the 
damage  it  produces,  and  I  shall  dis¬ 
cuss  briefly  how  destruction  occurs. 
Next  I  shall  describe  the  sequence  of 
events  that  leads  to  explosions  and 
the  processes  in  the  source  region 
during  the  explosion  itself.  Then  1 
shall  consider  the  current  status  of 
research  on  accidental  explosions. 

Damage  mechanisms 

Explosion  damage  occurs  for  a  num¬ 
ber  of  reasons.  The  enclosing  struc¬ 
ture  can  be  blown  apart  and  become 
the  source  of  primary  fragments.  The 
blast  wave  can  inflict  damage  on 
objects  it  encounters,  sometimes 
producing  secondary  fragments  that 
cun  damage  other  objects  (r.g.  flying 
glass  from  windows).  Radiation  from 
certain  combustion  explosions  can 
cause  burns  and  ignite  combustible 
material  at  a  distance. 


Internal  damage  to  the  enclosing 
structure,  though  difficult  to  quan¬ 
tify,  can  be  used  to  evaluate  the 
course  of  an  explosion.  In  buildings 
the  direction  in  which  objects  such  as 
walls,  partitions,  furniture,  etc.,  are 
displaced  can  usually  be  used  to 
diagnose  the  progress  of  the  explo¬ 
sion.  Metal  vessels  rupture  in  either 
a  brittle  or  ductile  mode,  depending 
on  the  particular  metal,  its  tempera¬ 
ture,  and  the  rate  of  application  of  the 
internal  pressure.  In  certain  cases, 
analysis  of  the  trajectories  of  primary 
fragments  can  be  used  to  deduce  the 
forces  acting  at  the  moment  of  rup¬ 
ture  (U.S.  Atomic  Energy  Commis¬ 
sion  196(5).  Baker  and  co-workers 
(1975, 1978)  have  shown  that  primary 
fragment  patterns,  though  somewhat 
irreproducible,  can  be  predicted  with 
reasonable  accuracy. 

The  blast  wave  produced  by  an  ex¬ 
plosion  is  called  a  free-field  blast  wave 
until  it  interacts  with  objects  in  its 
path.  Its  structure  is  quite  dependent 
upon  the  behavior  of  the  source  of  the 
explosion.  If  the  source  has  a  very 
high  energy  density  (units,  J  m"3) 
and  also  a  high  power  density  (units, 
J  m~3  s_1),  i.e.  releases  the  energy 
very  rapidly,  both  the  source  and  the 
blast  wave  are  said  to  be  ideal. 

Our  understanding  of  the  structure, 
scaling  laws,  and  damage  potential  of 
an  ideal  blast  wave  is  now  well  de¬ 
veloped  (Lee  et  al.  1969;  Baker  1973; 
Swisdak  1975).  Three  kinds  of  ex¬ 
plosions  that  yield  ideal  blast  waves 
are  point  source  explosions,  nuclear 
explosions,  and  condensed  phase 
detonations.  A  point  source  explosion 
is  a  mathematical  idealization:  a  fixed 
amount  of  energy  is  released  at  a 
point  in  an  infinite  uniform  atmo¬ 
sphere  in  an  infinitesimal  time.  Its 
energy  and  power  densities  art- 
therefore  infinite.  A  small  volume  of 
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Figure  1.  The  blast  wave  from  an  ideal  explo¬ 
sion  is  o  simple  pressure  wave  with  spherical 
symmetry,  the  leading  edge  of  which  Is  a  shock 
wave.  In  its  wake,  the  pressure  drops  to  a  sub- 
atmospheric  level  (negative  phase)  before  re¬ 
turning  to  the  ambient  level  due  to  recom¬ 
pression.  The  shape  of  the  wave  at  any  instant 
in  time  is  shown  at  the  top,  and  the  experience 
of  a  stationary  observer  os  the  blast  wave  passes 
is  shown  at  the  bottom.  (After  Baker  1973.) 
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air  in  the  vicinity  of  the  source  is 
compressed  and  heated,  thereby 
producing  a  region  of  very  high  pres¬ 
sure  that  is  bounded  by  a  strong 
spherical  shock  wave.  This  wave  de¬ 
cays  as  it  travels  away  from  the  source 
and  eventually  becomes  a  weak 
acoustic  wave. 

Nuclear  explosions  and  condensed 
phase  detonations  also  have  very  high 
energy  and  power  densities.  Their 
blast  waves  become  ideal  waves  be¬ 
fore  their  shock  pressures  have 
dropped  below  the  levels  required  for 
total  destruction  of  objects  in  their 
paths.  Thuu,  for  all  practical  pur¬ 
poses,  their  blast  waves  may  be  con¬ 
sidered  ideal,  even  though  there  are 
near-field  differences.  Sachs  (1944) 
was  the  first  to  show  that  the  tempo¬ 
ral  structure  of  all  ideal  blast  waves  is 
similar,  and  that  their  properties  are 
related  only  to  the  total  energy  of  the 
source  and  the  distance  of  the  ob¬ 
server  from  the  source. 

A  free-field  blast  wave  interacts  with 
objects  it  encounters  in  a  very  com¬ 
plex  manner.  The  shock  is  reflected 
and  refracted,  and  produces  higher 
pressures  locally  than  those  measured 
in  the  free  field,  thereby  causing 
objects  to  bend,  buckle,  tear,  or 
shatter. 


Figure  2.  For  lung  damage  in  man,  both  over¬ 
pressure  and  impulse  asymptotes  have  been 
observed.  For  eardrum  rupture,  only  the  ov¬ 
erpressure  asymptote  has  been  observed.  The 
response  of  the  eardrum  is  so  rapid  that  even 
a  blast  wave  with  a  small  impulse  appears  as  a 
pulse  of  long  duration,  making  it  impossible  to 
observe  the  impulse  asymptote.  K.ven  for 
temporary  hearing  impairment  (—24  hr),  it  has 
not  been  possihle  to  observe  the  impulse  as¬ 
ymptote,  even  though  the  curve  is  rising  toward 
the  asymptote.  (After  Raker  ct  al.  1973.) 


Sperrazza  (1951)  was  the  first  to  show 
empirically  that  blast  damage  in  any 
specific  mode  is  related  to  both  the 
peak  overpressure,  P,  and  the- positive 
impulse,  /,  of  the  wave  (Fig.  1).  Spe¬ 
cifically  he  observed  that  there  is  an 
impulse  asymptote  below  which  there 
is  no  damage,  irrespective  of  the 
magnitude  of  the  overpressure,  and 
an  overpressure  asymptote  below 
which  there  is  no  damage,  irrespective 
of  the  impulse  (or  duration)  of  the 
blast  loading.  An  example  of  this  be¬ 
havior  is  shown  in  Figure  2.  Newmark 
(1953)  independently  and  analytically 
identified  the  two  P-I  threshold 
asymptotes  for  damage,  using  the 
concept  of  energy  absorption.  Baker 
and  his  colleagues  (1978)  showed  re¬ 


cently  that  the  Jarrett  (1968)  formula 
for  blast  damage  to  brick  houses, 
which  was  developed  empirically 
from  bomb  damage  during  World 
War  II,  yields  typical  asymptotic  be¬ 
havior  when  plotted  in  the  P-I  plane. 
This  is  especially  interesting  since, 
unlike  Sperrazza,  they  did  not  spe¬ 
cifically  identify  impulse  as  being 
important. 

The  response  of  objects  to  drag  forces 
induced  by  the  blast  wave  is  more 
complex  than  that  predicted  by  sim¬ 
ple  P-I  response  curves  (Baker  et  al. 
1978).  This  is  because  the  kinetic 
energy  of  the  flow  associated  with  the 
blast  wave  imposes  a  drag  force  which 
can  cause  bending  of  thin  objects  such 
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as  flagpoles  or  light  standards,  or 
tumbling  or  gross  displacement  of 
unattached  bodies  such  as  trucks  or 
people  (Baker  et  al.  1973).  The  drag 
force  is  dependent  on  both  the  shape 
of  the  blast  wave  and  the  detailed 
shape  and  orientation  of  the  object. 
Penny  and  his  colleagues  (1970) 
based  their  determination  of  the 
TNT  equivalent  of  the  Hiroshima 
and  Nagasaki  atomic  bombs  on  such 
drag  forces. 

If  the  accidental  explosion  is  a  com¬ 
bustion  explosion  and  produces  a 
fireball,  radiation  damage  can  be  se- 


Typo  of  explosion 
Condensed  phase  detonations 


Combustion  explosions  of  gaseous  or  liquid 
fuels  in  enclosures 

Combustion  explosions  of  dusts  in  enclosures 

Bolling-liquid  expanding-vapor  explosions 
(BlEVEs) 

Unconflned  vapor-cloud  explosions 

Explosions  of  pressurized  vessels  containing 
nonreactlve  gaseous  materials 


and  may  eventually  provide  infor¬ 
mation  on  the  extent  of  the  damage 
for  specific  exposure  levels  and  rea¬ 
sonable  criteria  for  allowable  fireball 
exposure. 

Types  of  accidental 
explosions 

The  term  accidental  explosion  covers 
a  wide  spectrum,  and  obviously  no 
two  explosions  are  exactly  alike.  It  is 
nevertheless  possible  to  group  acci¬ 
dental  explosions  into  nine  major 
categories,  each  with  its  own  distinc¬ 
tive  characteristics  (Table  1). 


Typical  situations 

Chemical  reactors,  distillation  columns, 
separators,  factories  that  produce  high 
explosives  and  propellants 

Buildings,  ships,  tankers,  boilers, 
compressed-air  linas 

Coal  mines,  grain  elevators,  pharmaceutical 
industry 

Ductile  vessels  such  as  tank  cars  containing 
combustible  high  vapor-pressure  liquids 

Spills  of  highly  volatile  fuels,  chemical  plants 

Boilers,  pressurized  gas  tanks 


highly  sensitive  substance  has  acci¬ 
dentally  been  allowed  to  concen¬ 
trate. 

When  large  quantities  of  high  explo¬ 
sives  are  bundled  or  stored  in  bulk  or 
in  containers  that  are  in  relatively 
close  contact,  detonation  can  be  truly 
disastrous.  In  1921,  in  Oppau,  Ger¬ 
many,  a  congealed  mass  of  4.1  X  106 
kg  of  ammonium  nitrate-sulfate 
double  salt  that  was  being  broken  up 
with  dynamite  for  use  as  fertilizer 
detonated  as  a  unit,  killing  an  esti¬ 
mated  1,100  people,  causing  severe 
damage  up  to  6  km  away,  and  leaving 
a  crater  130  m  in  diameter  and  60  m 
deep.  Until  the  time  of  the  disaster,  it 
had  been  thought  that  this  salt,  even 
though  exothermic,  could  not  support 
a  detonation.  Even  today,  labora¬ 
tory-scale  testing  would  indicate  that 
it  is  nondetonatable  because  its 
minimum  charge  diameter  for  deto¬ 
nation  is  much  larger  than  the  size 
usually  used  for  sensitivity  tests. 
Another  major  explosion  involving 
ammonium  nitrate  took  place  in 
Texas  C'ty  in  April  1947,  when  the 
ship  Grand  Camp  caught  fire 
(Wheaton  1948).  The  fire  accelerated 
out  of  control,  causing  the  ammonium 
nitrate  in  the  hold  to  detonate.  All 
houses  within  a  1.5  km  radius  were 
totally  destroyed,  and  it  was  esti¬ 
mated  that  516  people  were  k;l!ed  and 
that  property  damage  amounted  to 
$67  million. 


Table  1.  Categories  of  accidental  explosions  and  where  they  are  most  likely  to  occur 


Explosions  resulting  from  chemical  reactor 
runaway 

Physical  vapor  explosions 


Explosions  resulting  from  nuclear  reactor 
runaway 


vere.  The  amount  of  radiation  is  pri¬ 
marily  a  function  of  the  size  and  du¬ 
ration  of  the  fireball,  since  fireball 
temperatures  for  most  fuels  are  ap¬ 
proximately  the  same  (»»  1350°K). 
Recently  it  has  been  postulated  that 
fireball  damage  can  be  evaluated  by 
a  technique  similar  to  the  P-I  tech¬ 
nique  for  blast-wave  damage  (Baker 
et  al.  1978).  In  this  case  the  static 
pressure  asymptote  is  replaced  by  a 
radiation  flux  rate  asymptote  for  long 
exposures,  and  the  impulse  asymp¬ 
tote  by  an  asymptote  representing  the 
total  flux  for  exposures  of  short  du¬ 
ration.  The  approach  looks  promising 
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Chemical  reactors 

Island  volcanoes,  molten  material  poured  into 
a  moist  container 

Nuclear  reactors 


Condensed  phase  detonations  ap¬ 
proximate  ideal  explosions — those 
with  infinite  energy  and  power  den¬ 
sities — although  in  practice  the  blast 
wave  is  usually  attenuated  somewhat 
due  to  confinement.  Blast  and  frag¬ 
ment  damage  can  be  estimated  if  the 
quantity  of  material  and  the  degree  of 
confinement  are  known  (AMCP  1972; 
Baker  1973;  Swisdak  1975).  There 
have  been  detonations  of  this  sort 
during  the  manufacture,  transport, 
storage,  and  use  of  high  explosives 
and  propellants,  and  in  chemical  re¬ 
actors,  distillation  columns,  separa¬ 
tors,  etc.,  when  some  unwanted  and 


The  second  type  of  explosion,  the 
combustion  explosion  of  a  gaseous  or 
liquid  fuel  in  an  enclosure,  has  two 
distinct  limit  behaviors.  An  enclosure 
that  has  a  length-to-riiameter  ratio 
(L/D)  of  about  one  and  that  is  not  too 
cluttered  with  equipment,  partitions, 
etc.,  will  usually  suffer  a  simple 
overpressure  explosion.  The  rise  in 
pressure  is  relatively  slow,  and  the 
weakest  windows  or  walls  blow  out 
first.  In  a  simple  frame  building,  the 
ceiling  will  rise  and  the  walls  will  all 
fall  out  at  about  the  same  time.  In  a 
steel  container  like  a  ship  hold  or  a 
boiler,  the  enclosure  will  tend  ;o  be¬ 
come  spherical  until  a  tear  or  rip 
vents  the  contents.  Although  the 
damage  to  the  enclosure  may  be  ex¬ 
tensive,  the  resultant  blust  wave  is 
ordinarily  quite  weak,  because,  in 
general,  buildings,  ships,  or  boilers 
vent  at  very  low  overpressures  (7-  70 
kPa  or  1-10  psig). 

In  enclosures  that  have  large  L/D  ra¬ 
tios  or  contain  large  pieces  of  equip- 


Figure  3,  An  internal  natural  gas  explosion  in 
the  elevator  shaft*  (located  on  the  left-hand 
face)  of  a  25-story  commercial  building  blew 
out  all  the  bricks  surrounding  the  shafts  and 
virtually  all  the  windows.  The  accident  oc¬ 
curred  in  April  1974  in  New  York  City.  (From 
NTSB  1976.) 


ment  or  internal  partitions,  flame 
propagation  following  ignition  causes 
gas  motion  ahead  of  the  flame,  which 
generates  turbulence  and  large-scale 
eddy  folding  where  the  flow  is  inter¬ 
rupted  by  obstacles.  This  in  turn 
causes  a  rapid  increase  in  the  effective 
flame  area,  resulting  in  a  more  rapid 
rise  of  pressure  and  further  turbu¬ 
lent/eddy  interactions.  Local  pres¬ 
sures  can  become  very  high  (*<1.5 
mPa  or  15  atm)  very  rapidly  (<1/ 
1,000  sec)  and  can  lead  to  gas-phase 
detonations  with  highly  localized 
massive  damage.  The  damage  is  usu¬ 
ally  greatest  furthest  away  from  the 
source  of  ignition.  These  detonations 
frequently  produce  strong  blast  waves 
and  high-velocity  fragments,  causing 
more  damage  to  the  surroundings 
than  simple  overpressure  explosions. 
Commercial  buildings  often  explode 
in  a  manner  indicative  of  significant 
wave  propagation  and  flame  acceler¬ 
ation  (Fig.  3).  Explosions  of  this  sort 
also  take  place  in  compressed-air 
lines,  where  the  fuel  is  oil  or  char  on 
the  walls  (Burgoyne  and  Craven 
1973). 

Halvorsen  (1975)  presents  examples 
of  tanker  explosions,  most  of  which 
are  simple  pressure  explosions  that 
occur  in  spaces  with  low  L/D  ratios. 
However,  some  tanker  explosions 
produce  the  highly  localized  damage 
typical  of  explosions  in  high  L/D  en¬ 
closures.  Figure  4  shows  the  remains 
of  such  an  explosion,  caused  by  ex¬ 
ternal  ignition  of  escaping  fuel  vapor 
during  ballasting.  Many  tanker  and 
supertanker  explosions  are  ignited  by 
sparks  of  static  electricity  generated 
by  the  high  pressure  water  spray  used 
for  cleaning. 

Combustion  explosions  of  dusts  in 
enclosures  can  be  quite  disastrous 
(Palmer  1973;  Bartknecht  1978). 
Contrary  to  some  commonly  held 
beliefs,  virtually  all  organic  dusts,  as 
well  as  certain  inorganic  or  metallic 
dusts,  uro  combustible  in  air  and  can 
explode  if  enclosed.  Dust  explosions 
exhibit  the  same  L/D  limit  behaviors 
as  gas  and  vapor  explosions,  but  the 
sequence  of  events  is  different. 


Combustible  quantities  of  airborne 
dust  normally  exist  only  inside  duct¬ 
ing  or  process  equipment  (e.g.  ham¬ 
mer  mills,  driers,  etc.).  A  small  ex¬ 
plosion  in  such  equipment  causes  it  to 
rupture  and  throw  burning  dust  into 
the  work  place.  If  the  work  area  is 
dirty,  the  resulting  gas  motion  and 
the  vibration  of  equipment  cause  the 
layered  dust  to  become  airborne  and 
then  fuel  a  more  violent  second  ex¬ 
plosion,  which  travels  through  the 
work  place  causing  major  damage. 

In  another  typical  sequence,  a  pile  of 
dust  starts  to  smolder  either  by 


spontaneous  combustion  or  because 
it  is  covering  a  hot  object  such  as  a 
motor  housing  or  a  lamp  fixture.  A 
worker  finds  the  fire  and  attempts  to 
put  it  out  with  a  chemical  extin¬ 
guisher  or  a  water  hose,  stirring  up  a 
large  cloud  of  dust,  a  portion  of  which 
is  already  burning,  and  an  explosion 
results. 

Coal  mine  explosions  have  been 
common  since  the  start  of  the  indus  ¬ 
trial  revolution.  Most  industrial  na¬ 
tions  have  extensive  research  pro¬ 
grams  to  investigate  such  explosions, 
and  the  dynamics  of  the  explosion 
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process  are  reasonably  well  under¬ 
stood.  Ignition  is  usually  effected  by 
the  release  and  ignition  of  a  pocket  of 
methane  in  the  mine,  or  by  detona¬ 
tion  of  the  explosives  used  in  blasting. 
Since  a  mine  is  always  quite  dirty, 
because  it  is  expensive  to  remove  all 
the  coal  dust,  the  resulting  dust  ex¬ 
plosion  propagates  over  long  dis¬ 
tances.  Although  ihe  annual  number 
of  deaths  from  coal  mine  explosions 
in  the  United  States  has  decreased  by 
more  than  half  since  the  beginning  of 
the  century  (Cyhulski  1975),  mine 
explosions  will  undoubtedly  occur  in 
the  future. 


Figure  4,  While  the  t.iherian  tnnker  .S\.S\ 
Sansiru'na  (250  m  long  and  52  m  wide)  was 
being  ballasted  in  I.os  Angeles  harbor  in  De¬ 
cember  1976,  a  cloud  of  combustible  vapor 
from  the  hold  formed  on  the  deck  and  was  ig¬ 
nited  by  an  unknown  source.  The  resulting 
flash  fire  entered  the  hold,  setting  of!  an  ex¬ 
plosion.  Witnesses  reported  that  the  entire 
tank  deck  and  midship  deckhouse  rose  250  m 
into  the  air,  leaving  only  the  part  shown  here. 
Severe  damage  was  reported  up  to  2  km  away, 
and  windows  were  broken  up  to  4  km  away. 
The  explosion  killed  6,  left  3  missing,  injured 
58,  und  caused  $21.6  million  in  damage.  (From 
IJSCC.  1977.) 


Dust  explosions  iti  grain  elevators  and 
in  various  milling  industries  are  also 
common,  averaging  ten  to  fifteen  a 
year  in  the  United  States.  In  Decem¬ 
ber  1977  two  disastrous  explosions 
took  place  within  five  days  of  each 
other,  at  the  Continental  grain  facility 
in  Westwego,  DA.  and  at  the  Farmers 
Export  grain  elevator  in  Galveston, 
TX,  significantly  reducing  the 
grain-exporting  capacity  of  the 
United  States. 

Boiling-liquid  expanding-vapor  ex¬ 
plosion x,  or  BLEVEs,  occur  when  a 
ductile  vessel  containing  a  liquid  with 


a  vapor  pressure  well  above  atmo¬ 
spheric  pressure  tears  open.  In  these 
explosions  the  tearing  process  is  rel¬ 
atively  slow,  and  a  small  number  of 
large  fragments  are  produced.  Be¬ 
cause  they  are  backed  by  a  liquid  that 
can  evaporate  very  rapidly  (flash 
evaporate),  even  large  fragments  can 
be  thrown  up  to  700  in.  If  the  liquid  in 
the  tank  is  combustible  and  the 
BLEVE  is  caused  by  heat  from  an 
external  fire,  a  buoyant  fireball  is 
produced,  the  duration  and  size  of 
which  are  determined  by  the  total 
weight  of  fluid  enclosed  in  the  tank. 
The  BLEVE  of  a  standard  311,000 
gallon  tank  car  containing  a  flash¬ 
evaporating  fuel  such  as  liquid  pe¬ 
troleum  gas  will  create  a  fireball,  the 
radius  of  which  is  ~50  m  at  ground 
level  and  ~100  m  once  it  becomes 
airborne.  Thus,  even  when  the  blast 
damage  is  mild,  missile  and  fireball 
damage  can  he  catastrophic. 

In  an  unconfined  vapor-cloud  ex¬ 
plosion,  the  first  event  is  the  massive 
spill  of  a  combustible  hydrocarbon 
into  the  atmosphere.  Following  this, 
one  of  four  things  can  happen:  the 
spill  may  he  dissipated  harmlessly 
without  ignition;  it  may  be  ignited 
immediately  upon  release,  causing  a 
fire  hut  no  explosion;  it  may  he  dis¬ 
persed  over  a  wide  area  before  being 
ignited,  causing  a  large  fire;  or,  if  the 
flame  accelerates  rapidly  enough,  a 
dangerous  blast  wave  can  he  pro¬ 
duced  (Strehlow  1973). 

Urcunfined  vapor-cloud  explosions 
can  he  both  very  spectacular  and  very 
dangerous,  because  the  leak  is  into 
the  open  air,  and,  with  the  right  me¬ 
teorological  conditions,  extremely 
large  clouds  of  the  combustible  mix¬ 
ture  can  he  produced  before  ignition 
occurs.  In  dune  1974  at  the  Nypro 
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chemioal  plant  near  Flixborough, 
England,  approximately  46,009  kg  of 
cyclohexanone  at  a  pressure  of  850 
kPa  (125  psig)  and  a  temperature  of 
155°C  was  released  through  two  0.7  m 
diameter  stub  pipes  and  flash-evap¬ 
orated,  producing  a  sizable  cloud 
throughout  the  plant.  Ignition  prob¬ 
ably  occurred  at  a  furnace  in  the  hy¬ 
drogen  plant  some  distance  from  the 
release  point,  and  the  fire  was  already 
large  by  the  time  the  flame  had  ac¬ 
celerated  enough  to  produce  a  blast 
wave  that  extensively  damaged  the 
plant  and  houses  up  to  2  km  away 
(Gugan  1978). 

An  incident  in  Franklin  county,  MO, 
in  1970  is  unique  in  that  the  vapor 
cloud  appears  to  have  detonated  as  a 
unit.  Because  a  considerable  amount 
of  the  fuel-air  mixture  was  too  rich  to 
detonate,  a  fireball  ensued,  burning 
an  area  of  14  ha. 

Explosions  of  pressurized  vessels 
containing  nonreactive  gaseous  ma¬ 
terials  cause  massive  damage  if  the 
contents  are  at  a  very  high  pressure. 
Although  the  fragments  are  not 
thrown  nearly  as  far  as  they  are  with 
a  BLEVE,  the  blast  is  more  severe, 
because  the  entire  contents  of  the 
vessel  contribute  to  the  blast  wave.  If 
the  vessel  is  weak,  it  will  rupture  at 
relatively  low  pressures.  A  large-scale 
utility  boiler,  for  example,  will  “ex¬ 
plode”  when  the  internal  pressure 
rises  only  1-15  kPa  (1.5-2  psig).  This 
pressure  rise  can  be  triggered  by  a 
combustion  explosion  in  the  boiler  or 
by  a  massive  steam  leak  caused  by  the 
rupture  of  a  large  tube  or  header.  In 
this  case  steam  enters  the  red  hot 
boiler  so  fast  that  the  normal  open¬ 
ings  are  insufficient  to  stop  the  pres¬ 
sure  rise.  In  either  case,  because  of  the 
low  L/D  ratio  of  the  boiler  and  its 
weak  structure,  the  damage  to  the 
surroundings  is  slight,  even  though 
the  boiler  itself  may  be  severely 
damaged. 

Chemiccl  reactor  runaway  may  re¬ 
sult  from  an  imbalance  in  the  con¬ 
trolled  reaction  through  loss  of  ade¬ 
quate  cooling,  inadequate  stirring,  the 
presence  of  too  much  catalyst,  etc. 
The  pressure  buildup  is  considerably 
slower  than  it  would  be  in  an  explo¬ 
sion  due  to  the  detonation  of  the 
contents  of  the  container,  and  the 
vessel  usually  ruptures  in  a  ductile 
mode.  If  the  contents  are  gaseous,  the 
rupture  resembles  that  of  an  over¬ 
pressurized  vessel.  If  the  contents  are 


liquid  and  above  the  flash-evapora¬ 
tion  temperature,  as  is  commonly  the 
case,  the  explosion  resembles  a 
BLEVE. 

Physical  vapor  explosions  occur 
when  a  finely  divided,  hot,  solid  ma¬ 
terial  is  rapidly  mixed  with  a  much 
cooler  liquid,  or  when  two  liquids  at 
different  temperatures  are  violently 
mixed  (Reid  1976).  The  cooler  liquid 
is  converted  to  vapor  sc  fast  that  an 
explosion  results  from  the  localized 
high  pressures  that  are  produced. 
These  explosions  have  been  observed 
when  molten  material  is  poured  into 
a  container  that  is  moist,  and  when 
cryogenic  liquid  natural  gas  that 
contains  approximately  10%  of  a 
higher  hydrocarbon  is  spilled  on 
water.  There  is  some  indication  that 
the  mixing  of  seawater  with  hot 
magma  may  have  been  responsible  for 
the  catastrophic  explosions  of  island 
volcanoes  such  as  Krakatoa  in  1883 
and  Surtsey  in  1963. 

Many  scenarios  have  been  con¬ 
structed  for  the  consequences  of  nu¬ 
clear  reactor  runaway  or  core  melt¬ 
down,  including  catastrophic 
breaching  of  the  containment  vessel 
by  an  internal  combustion  explosion 
or  a  simple  pressure  burst.  Nuclear 
reactors  are  constructed  in  such  a  way 
that  an  accidental  upset  could  not 
possibly  produce  anything  remotely 
resembling  an  atomic  bomb  explo¬ 
sion.  If  a  nuclear  reactor  were  to  ex¬ 
plode,  the  attendant  release  of  long- 
lived  radioactive  material  would  be  so 
disruptive  to  the  local  area  that  the 
explosion  damage,  no  matter  how 
extensive,  would  be  considered  minor. 
Fortunately,  no  such  accident  can  be 
cited  because  none  has  occurred,  an 
enviable  and  unequaled  safety  record 
for  the  first  twenty  years  of'  mer- 
cialization  of  a  new  technology. 

Research  on  accidental 
explosions 

Research  on  accidental  explosions  has 
been  directed  toward  determining 
how  the  blast  waves  they  produce 
differ  from  ideal  blast  waves,  eluci¬ 
dating  the  explosion  process,  and 
devising,  developing,  and  proofing 
techniques  to  mitigate  their  efiects 
and  reduce  their  frequency.  Only  re¬ 
cently  have  nonideal  explosions  been 
studied  in  any  detail  (Strehlo\<  and 
Baker  1976).  It  has  been  shown  that 
nonideal  behavior  of  a  source  is  al¬ 
ways  directly  responsible  for  any 


nonideal  behavior  of  the  resulting 
bla9t  wave.  A  source  is  nonideftl  if  it 
has  a  low  energy  or  power  density,  if 
the  release  of  energy  is  nonspherical, 
or  if  it  is  totally  or  partially  confined 
by  inert  material. 

Explosions  of  pressurized  vessels  and 
combustion  explosions  of  gases,  va¬ 
pors,  and  dusts,  in  or  out  of  enclo¬ 
sures,  are  examples  of  explosions  of 
low  energy-density  sources.  The  ef¬ 
fect  of  low  energy  density  on  the 
structure  of  blast  waves  has  been 
studied  rather  extensively  for  burst¬ 
ing  spheres,  both  experimentally  and 
theoretically.  Figure  5  shows  the  ev¬ 
olution  of  the  blast  wave  from  such  a 
bursting  sphere. 

The  effect  of  slow  energy  addition, 
and  the  concomitant  low  power  den¬ 
sity,  have  been  studied  for  two  cases, 
also  illustrated  in  Figure  5:  Adamczyk 
(1976)  studied  ramp  addition  of  en¬ 
ergy  (the  addition  of  energy  at  a  uni¬ 
form  rate  throughout  the  source  re¬ 
gion),  whereas  Kuhl  and  Strehlow 
and  their  co-workers  (Kuhl  et  al. 
1973;  Strehlow  et  al.  1979)  studied 
flame  addition  of  energy  (by  central 
ignition  in  a  spherical  source  region 
initially  at  atmospheric  pressure).  In 
a  very  early  piece  of  work,  Taylor 
(1946)  predicted  analytically  the  flow 
structure  for  a  spherical  piston  ex¬ 
panding  at  constant  velocity. 

It  has  been  found  that  for  spherical 
source  regions  low  energy  and  power 
densities  generally  cause  the  over¬ 
pressure  in  the  blast  wave  to  fall 
below  that  predicted  by  Sachs’s 
scaling  on  the  basis  of  the  energy  of 
the  source.  It  has  also  been  shown 
that  all  such  nonideal  sources  yield 
blast  waves,  the  positive  phase  im¬ 
pulse  of  which  is  equal  to  that  of  an 
ideal  explosion  of  the  same  total  en¬ 
ergy. 

Nonspherical  or  directional  blast 
waves  are  generated  by  the  firing  of  a 
gun  or  rocket,  the  cracking  of  a  whip, 
and  lightning.  There  have  been  a  few 
limited  analytic  studies,  a  few  nu¬ 
merical  studies,  and  some  very  spe¬ 
cific  experimental  st  udies  on  the  ef¬ 
fect  of  nonspherical  release  of  energy 
on  the  structure  of  the  blast  wave 
produced.  Unfortunately,  the  analytic 
studies  have  all  been  based  on  very 
specialized  geometry  or  the  use  of 
approximations  that  tend  to  be 
physically  untenable,  while  the  nu¬ 
merical  studies  have  required  two- 
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and  three*  dimensional  hydrocodes 
that  are  extremely  expensive  to  run. 
Thus,  generalizations  are  not  cur¬ 
rently  available.  1  recently  discovered 
that  the  overpressure  generated  by 
the  deflagration  of  a  cloud  of  arbi¬ 
trary  shape  can  be  estimated  by 
applying  a  fundamental  acoustic 
principle  first  stated  by  Stokes  in 
1849  (Strehlow,  unpubl.),  but  the  full 
implications  of  this  observation  have 
yet  to  be  examined. 

Confinement  weakens  the  external 
blast  wave  and  causes  it  to  become 
highly  directional.  Also,  a  confined 
explosion  may  cause  massive  internal 
damage  but  relatively  little  external 
damage,  as  in  explosions  in  coal  mines 
and  some  buildings. 

Overpressures  close  to  the  source  are 
much  lower  for  nonideal  than  for 
ideal  explosions.  If  an  explosion  is 
deflagrative,  no  shock  wave  is  formed 
in  the  near  field  (Pig.  5).  It  was  first 
noted  by  Rayleigh  (1876)  that  explo¬ 
sions  with  low  power  densities  pro¬ 
duce  a  blast  wave  with  a  sizable  neg¬ 
ative  phase  impulse.  He  inferred  this 
because  he  sometimes  found  the  glass 
from  a  shattered  window  of  a  nearby 
building  outside  the  building. 


The  explosion  process  is  reasonably 
well  understood  except  for  uncon¬ 
fined  vapor-cloud  explosions.  What 
is  not  yet  understood  is  how,  in  cases 
of  soft  ignition  (i.e.  an  open  flame 
such  as  a  pilot  light),  an  essentially 
unconfined  flame  accelerates  suffi¬ 
ciently  to  produce  a  damaging  blast 
wave.  Prior  to  1975,  transition  to 
detonation  was  thought  to  occur  only 
if  a  portion  of  the  gas  was  heated 
rapidly  by  compression  to  its  auto¬ 
ignition  temperature.  Then  in  1976, 
an  accidental  detonation  took  place 
in  Wagner’s  laboratory  in  Gottingen, 
under  conditions  in  which  compres¬ 
sion  heating  was  impossible  (Wagner, 
pers.  comm.).  Sinco  then,  Knystautas 
and  co- workers  (1979)  have  shown 


that  shockless  initiation  is  possible  in 
large-scale  eddies  formed  near  ob¬ 
stacles  in  the  flame-induced  flow 
ahead  of  the  flame. 

Research  into  the  detonability  of 
natural  gas-air  mixtures  is  important, 
because  liquid  natural  gas  is  being 
shipped  in  bulk  (25,000  m3  tanks,  5  to 
a  tanker),  and  accidents  have  the 
potential  for  producing  extremely 
large  vapor  clouds.  It  is  known  that, 
direct  initiation  of  detonation  in  an 
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unconfined  cloud  of  pure  methane  collection  systems)  is  much  more 
requires  a  very  large  charge  of  a  high  difficult.  It  is  unnecessary  with  gases 
explosive,  but  how  large  is  still  a  and  vapors  if  enough  air  is  injected  to 
matter  of  contention.  The  problem  is  dilute  the  fuel  to  a  concentration  of 
complicated  by  the  fact  that  natural  less  than  one-quarter  of  its  lean 
gas  usually  contains  a  small  percent-  flammability  limit  (i.e.  the  largest 
age  of  higher  hydrocarbons,  and  the  concentration  of  fuel  and  air  at  which 
addition  of  about  10%  ethane  to  a  flame  will  just  fail  to  propagate  over 
methane  is  known  to  increase  its  long  distances,  which  is  about  1%  fuel 
sensitivity  to  direct  initiation  con-  in  air  for  most  higher  hydrocarbons), 
siderably.  Dusts  present  a  special  problem,  be¬ 

cause  if  the  flow  velocity  at  any  point 
Venting  has  been  found  to  be  a  viable  in  the  system  falls  below  a  critical 
way  of  limiting  damage  from  explo-  value,  the  dust  settles  in  the  duct, 
sions  in  enclosures  containing  com-  thereby  producing  a  dangerous  sit- 
bustible  gases,  vapors,  or  dusts,  and  uation.  Requirements  for  venting 
from  runaway  chemical  reactors  ducts  are  given  by  the  National  Fire 
(Singh  1978).  Bradley  and  Mitcheson  Protection  Association  (1978).  The 
(1978a, b)  recently  constructed  a  principle  of  chemical  reactor  venting 
rather  complete  theory  of  venting,  using  a  burst  disk  is  well  understood 
and  their  theoretical  estimates  ap-  if  the  reactor  contains  a  gas,  but  not 
proximate  closely  virtually  all  the  if  it  contains  a  flash-evaporating  liq- 
extant  laboratory  results.  However,  it  uid,  because  of  the  difficulty  of  cal- 
has  been  shown  that  for  gas-air  ex-  culating  the  mass  flow  and  the  rate  at 
plosions  in  room-size  vessels,  the  which  the  pressure  drops. 

Bradley  theory  underestimates  the 

maximum  internal  pressure.  Evi-  In  contrast  to  venting,  which  is  a 
dently  an  acoustic  instability  gener-  passive  technique  for  the  mitigation 
ated  by  flame  propagation  causes  of  the  effects  of  explosions,  suppres- 
higher  overpressures  than  predicted,  sion  can  be  either  passive  or  active. 
Such  instabilities  do  not  appear  in  The  water  barrier  method  used  to 
small  chambers  because  their  high  suppress  coal  mine  explosions  is  a 
natural  frequencies  preclude  direct  typical  passive  technique  (Liebman 
interaction  between  the  acoustic  and  Richmond  1974).  Most  coal  mine 
modes  and  the  flame.  explosions  originate  at  the  working 

face  of  the  mine  and,  as  is  chnracter- 
Dust  explosion  venting  has  been  istic  of  explosions  in  enclosures  with 
studied  extensively  by  Bartknecht  large  L/D  ratios,  do  not  become  really 
(1978).  He  finds  that  for  vessels  with  dangerous  until  they  have  traveled 
small  L/D  ratios  and  volumes  larger  some  distance  from  their  point  of  or- 
than  0.02  ma,  one  can  define  a  con-  igin.  The  water  barrier  consists  of  a 
stant  for  any  substance  which  is  pro-  number  of  large  open  tubs  of  water, 
portional  to  the  maximum  rate  of  which  are  suspended  near  the  ceiling 
pressure  rise  and  which  can  be  used  to  of  the  mine  a  short  distance  from  the 
determine  the  proper  vent  area.  This  working  face.  The  air  flow  generated 
approach  to  the  rating  of  the  explos-  by  the  incipient  explosion  upsets  the 
iveness  of  dust  has  recently  been  tubs  and  disperses  the  water.  This 
adopted  by  the  National  Fire  Pro-  airborne  water  acts  as  a  heat  sink  and 
tection  Association.  Bartknecht’s  quenches  the  flame.  Another  passive 
work  shows  definitively  that  the  suppression  technique  consists  of 
Hartman  bomb,  which  was  used  ex-  sprinkling  the  mine  with  rock  dust, 
tensively  for  testing  the  explosiveness  usually  calcium  carbonate  from 
of  dust  in  the  1960s  in  the  United  crushed  shells.  This  inert  powder  is 
States  (Dorsett  et  al.  1960),  is  too  picked  up  with  the  coal  dust  and 
small  to  evaluate  properly  the  rate  at  lowers  the  temperature  enough  to 
which  the  pressure  rises,  particularly  quench  the  flame  (Richmond  ct  al. 
for  dusts  with  a  very  rapid  pressure  1979). 
buildup.  Venting  dust  explosions  is 

quite  different  from  venting  gaseous  Active  suppression  devices,  produced 
explosions,  since  there  is  considerably  commercially  in  the  United  States 
more  combustion  outside  the  enclo  and  West  Germany,  make  use  of 
sure.  pressurized  bottles  containing  a  dry 

powder  or  a  halogenated-hydrocar- 
In  general,  venting  structures  with  bon  extinguishing  agent,  which  are 
large  L/D  ratios  (such  as  the  ducting  triggered  to  open  by  a  sensor  that  is 
for  spray-painting  hoods  or  for  dust-  sensitive  to  either  a  rise  in  pressure  or 


.  the  light  cn  titted  by  a  flame  (Liebman 
•tat.  1979) 

One  way  of  ensuring  that  an  enclosure 
wilt  not  explode  is  to  make  it  inert 
with  nitrogen  or  some  other  gas.  After 
they  are  dried,  combustion  products 
from  boilers  are  both  a  suitable  and  a 
relatively  inexpensive  source.  The  one 
danger  always  present  is  the  possi¬ 
bility  of  suffocation  if  someone  enters 
the  enclosure. 

The  incidence  of  BLEVEs  can  be  re¬ 
duced  by  the  application  of  intum- 
escent  coatings  (Cagliostro  et  al. 
1975).  Such  coatings  swell  and  char 
readily  when  heated  by  an  external 
fire,  thereby  reducing  heat  transfer  to 
the  substrate,  extending  the  fire  life 
of  the  object,  and  reducing  the  rate  at 
which  heat  is  transferred  to  the  inte¬ 
rior.  Although  there  are  many  ade¬ 
quate  intumescent  formulations,  very 
few  are  rugged  enough  to  withstand 
normal  operations  of  a  tank  car  and 
the  shock  of  derailment  without  sep¬ 
arating  from  the  surface  and  leaving 
large  areas  unprotected.  However, 
some  promising  new  intumescent 


materials  are  now  being  tested  on 
tank  cars  in  regular  service  (NTSB 
1979). 

A  study  made  some  years  ago  showed 
that  virtually  all  tank  car  punctures 
resulting  from  accidents  took  place  in 
the  bottom  third  of  one  end  of  the  car 
and  were  caused  principally  by  the 
coupler  on  the  next  car.  In  1974  the 
United  States  government  first  re¬ 
quired  that  both  ends  of  tank  cars 
carrying  certain  hasardous  materials 
be  fitted  with  nonslip  couplers  and 
heavy  vertical  steel  plates  to  protect 
them  from  such  punctures.  A  recent 
study  (NTSB  1979)  shows  that  this 
rather  inexpensive  retrofit  program 
is  proceeding  according  to  schedule, 
the  only  major  problem  being  the 
slow  response  of  some  regulatory 
agencies. 

Current  status  and  future 
needs 

Our  understanding  of  accidental  ex¬ 
plosions  has  increased  tremendously 
in  the  past  fifteen  years,  primarily  as 
a  result  of  societal  concern.  In  the 


United  States,  for  example,  the  last 
decade  or  so  has  seen  new  and  im¬ 
portant  involvement  in  explosion 
protection  by  the  federal  government, 
insurance  underwriters,  trade  asso¬ 
ciations,  and  technical  societies.  De¬ 
spite  our  increased  understanding, 
however,  we  have  not  as  yet 
reached— and  never  will  reach— a 
safety  level  at  which  we  are  entirely 
free  of  risk  from  accidental  explo¬ 
sions.  If  we  as  a  society  are  dedicated 
to  reducing  the  risk,  we  must  continue 
to  support  research  on  the  causes  and 
mechanisms  of  explosions  as  well  as 
on  protective  devices. 

There  is,  moreover,  another  aspect  to 
this  problem  that  is  equally  impor¬ 
tant.  It  concerns  the  level  of  training, 
motivation,  and  general  under¬ 
standing  of  the  people  who  live  and 
work  in  our  highly  technological  so¬ 
ciety.  At  least  60%  of  all  accidental 
explosions  are  caused  by  human  error 
and  not  by  just  mechanical  malfunc¬ 
tion  of  equipment.  How  do  you  mo¬ 
tivate  a  chemical  plant  worker  to  re¬ 
frain  from  smoking  a  cigarette  behind 
a  large  oil  tank  where  his  fellow 


Figure  5.  The  spatial  distribution  of  pressure  at  equal  time  intervals  is 
shown  for  the  blast  waves  produced  by  three  nonideal  sources.  The  center 
of  the  source  is  to  the  right  in  each  case.  For  a  bursting  sphere  with  an 
initial  pressure  ratio  of  9  (/e/O,  there  is  immediate  formation  of  a  shock 
wave  and  later  formation  of  a  high-pressure  pulse  at  the  center  and  a 
resultant  second,  weaker  shock.  (After  l.uckritz  1977.)  The  ramp  addition 
of  energy  ( upper  right )  causes  a  compression  wave,  which,  given  suffi¬ 
ciently  high  energy  and  power  densities,  slowly  steepens  to  a  shock  wave 
as  it  propagates  away  from  the  source.  In  this  case  there  is  no  pressure 
pulse  at  the  center.  (From  Adamczyk  1976.)  For  a  flame  propagating  away 
from  the  center  of  a  spherical  source  region  at  constant  velocity  (/on  er 
right),  there  is  no  lt>ud  shock  wave  in  the  How  field  for  flame  velocities 
less  than  —35  m/sec.  The  oscillations  in  the  central  region  are  damped 
numerical  instabilities  due  to  the  low  rate  of  energy  addition  and  are  not. 
physically  real.  (From  Strehlow  et  al.  1979.) 


» 


worker?  cannot  see  him?  How  do  you 
train  workers  to  respond  appro¬ 
priately  in  an  emergency  situation? 
How  do  you  teach  all  the  people  who 
sweep  down  floors  and  walls  in  grain 
elevators  that  grain  dust  suspended 
in  air  is  dangerously  explosive?  How 
do  you  get  people  to  write  regulations 
that  reduce  risk  without  imposing  an 
excessive  economic  burden?  These 
questions  and  many  more  like  them 
must  be  addressed. 

It  appears  that  the  awareness  devel¬ 
oped  over  the  past  fifteen  years  will 
continue  to  grow.  I  can  see  no  dimi¬ 
nution  in  the  research  effort,  and  I  am 
heartened  by  the  large  number  of 
companies  and  government  research 
agencies  that  are  worried  about  safety 
in  the  work  place  and  in  transporta¬ 
tion  systems.  Accidental  explosions 
will  certainly  continue  to  occur,  and 
only  a  continued  effort  to  further  our 
understanding  will  lead  to  a  reduction 
in  their  frequency  and  impact  in  the 
future. 
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appendix  c 


Blast  Wave  from  Deflag  rati  ve 
Explosions:  an  Accoustic 
Approach 


r 

‘i  . 


Simple  acoustic  source  theory  has  been  applied  to 
determine  the  maximum  overpressure  Obtainable  by 
the  deflagration  of  nonsphericaj  Giouds,!Jn  three 'di¬ 
mensions,  Over-pressure  is  generated  not  by  the  rate  of 
energy  addition  but  by  the  first-tune  derivative  of  the 
‘rate.  'Therefore,  defiagrative  combustion  of  edge- 
ignited  clouds  produces  markedly  less  overpressure 
than  central,  spherical  ignition.  ' 

''  R.A.  Strehlow  ;  . 

University  of  Illinois  at  Urbana-champaigh,  Urbana,  IL. 


The  determination  of  the  blast  wave  produced 
in  free  space  by  the  nonspherical  combustion 
of  a  nonspherical  vapor  cloud  has  not. been 
examined  in  any  detail  .  either  analytically, 
or  numerically.  The  problem  is  important 
because  at  the  present  time  spherica'  tneory 
is  being  used  to  estimate  the  overpressure, 
for  clouds  of  arbitrary  shape  and  no  one 
knows  how  conservative  this  approach  is. 

This  paper  applies  simple  acoustic  source 

theory  tc  the  deflagrative  combustion  of  a 

few  simple  free,  cloud  shapes  to  determine  & 
first  cut  estimate  of  the  maximum  overores^ure 
that  one  could  expect  as  a  function  of  cloud 
site  and  shape. 

THE  SIMPLE  ACOUSTIC  SOURCE 

n ) 

Stokes  first  showed  that  a  simple 
source  of  mass  at  a  point  <h(t)  (Units:  Kg/ sec) 
generates  a  sound  wave  in  three  dimensions 
which  has  an  overpressure  which  is  proportion¬ 
al  to  ft(t).  Specifically,  as  Lighthill^! 
states, 


P  -  ? 


m 


t  - 


/4;rr 


(1) 


where  t  -  3—  replaces  t  because  the  wave  is 
oropagating°away  from  the  source  region  at 
the  velocity  of  sound,  a0.  Lighthiil  also 
states  that  at  some  distance  from  the  source 
toe  sound  will  appear  to  emanate  from  a  ooint 

-1^.S701-.8V3810  $02.00  !  ’981  MC,n£ 


(1  .a,.-,  'the  wave  will  be  very  close  to  .spheri¬ 
cal)  if  the  quantity  (wr/a0)2  <<1.  Here,  u> 
is  the  circular, frequency  of  the  source  of 
sPund  and  r  is  a  characteristic  radius  uf  . 
the  source  region.  Finally,  for  a  source 
which  has  a  finite  duration  Jfii{t)dt  =0, 
and  therefore  /(P-Po)dt  -  0  -in.  the  wave. 

This  means  that  the  acoustic  pulse  from  a 
source  of  finite  duration  must  have  a' 
negative  phase  impulse  equal  to  the  positive 
phase,  impulse.'  This  should  b.e  .contrasted 
to  the  pulse7' emitted  by  a  finite  duration 
source  in  a  strictly  one-dimensional  channel.’ 
In  this  case  the  acoustic  overpressure  :is 
proportional  to  ih(t)  and  the  overpressure  ,p 
in  the  pulse  is  always  positive'  everywhere.'  ■* 

how  consider  the  deflagrative  combustion 
of  an  unconfined  cloud  of  arbitrary  shape. 

For  simple  source  acoustic  behavior  the 
effective  rate. of  mass  addition  fi(t)  may  be 
replaced  by.  an  effective  rate  of  volume 
addition  multiplied'  by  the  initial  gas 
density,  that  is,  m(t)  *  o0$(t). 


However,  for  deflagrative  combustion 

V.  -  v  ^ 

m(t)  *  PQv(t: .  “ ‘ 


4r(S  (ti  ■  A.(t) ) 


where  S u ( t }  is  tne  effective  normal  burning 
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velocity.,  Af(t)  is  the  effective  frontal  area 
of  the  flame  and. the  ratio  C Vb-Vy ) /V u _is  the 
new  volume  produced  per  unit  volume  of  gas 


burned  at  constant  pressure.  Using  the 
notation  of  Strehlow  et  al .  ill  (Vb-'fj i/'fj  s 


an  effective  dimensionless 


q/vA.  Here,  A  is  u  ,  . 

neat  addition.  Its  relation  to  the  o,  or 
.Strenlow  atal .  is  given  in  Appendix  A. 
Therefore,  using  the  definition  of  the  vel¬ 
ocity  of  ■  sound  in  an  ideal  gas  and  combining 
Eqs.  ,(!■)  and  (2)  yiekis 


-  v  'ft'"-  d. 

0  S'-  -  —  - : — 


*  CSJF  It  lsu(t> 


AfU)J 


the  relation 


*  is,  dAf 


9  = 


A-a^r  dt 
o 


(A) 


vrere  ^  is  the  distance  frori  the  source_t‘ 
‘■-•e  observer..  For  a  growing  s chemical  * •; 
Sal  1 


Af  ■  4tt£ 


and 


6r<  -  s  .  s  -i. 
S 


dt 


Substituting  yields 


where  P  is  the  acoustic  overpressure. 

Equation  (3)  is  a  significant  relationship. 

It  states  that  the  acoustic  overpressure 
1  genfe'ritiad'vqy.  deflagrative  combustion  in 
three  dimensions  is  given  by  the  time  rate 
of  change  of  the  rate  of  energy  addition 
rather  than  the  rate  of  energy  addition  as 
pit  is  in.  strictly  one-0  flow. 

The  deflagrative  combustion  of  a  cloud 
generates  a  sound  wave  of  very  low  frequency. 
The  period  of  the  sound  source  in  this  case 
is  given  by  x  =  2m/ui  where  t  is  the  burn 
time...  Thus  the  requirement  that  the  source 
region  is  acoustically  simple  becomes 


P  = 


) 


ofr 


V 


where  the  r  in  the  denominator  must  be  equal 
to  or  greater  than  rf.  If  we  let  it  be  equal 
to  rf,  and  if  we  substitute  the  definition 
of  q  given  above,  we  obtain  ah  expression 
for  the  acoustic  pressure  on  the  unburned 
side  of  the  flame 


r  \  t 


[2ir(rm/T)/a0]2  «  1 

where  rm  is  a  major  dimension  of  the  cloud 
and  approximates  the  distance  traversed  by 
the  deflagration  wave.  Even  with  distortion 
and  motion  due  to  flame  propagation  r/x  ~  Su 
and  the  simple  source  requirement  becomes 
4xt2M|  <<  1.  If  the  quantity  is  to  be  less 
than  Q.l  for  no  blast  wave  distortion  Su 
‘.'win  be  about  17.5  m/sec  before  deflagrative 
behavior  produces  significant  distortion  of 
the  blast  wave.  Even  for  significantly 
higher  values  of  Su  the  first  effect  will  be 
only  to  distort  the  shape  of  the  blast  wave. 
Higher  overpressures  than  estimated  here  will 
not  occur  until  Sm  becomes  quite  large. 

This  limitation  or  the  theory  has  not  been 
quantified  as  yet. 

APPLICATION  IN  SPHERICAL  GEOMETRY 


which  is  identical  to  Eq.  (19)  of  Strehlow 
et  i!_*  *  ~  ■  ".V 

If  we  write  rf  =  S^t  for  the  flame 
propagation  in  Eq.  (5)  and  recall  that  away 
from  the  flame  t  becomes 


Eq.  (5)  may  be  written 


(7) 


P  =  2v 


1  - 


V  ' 3  ' 

_£  3 

V  i  'Su 


Uj 


at  1 

-f -’j 


However,  the  transformation  from  a 
flame  to  a  Taylor  spherical  oiston  is 


Two  classical  problems  have  been  treated 
in  spherical  geometry.  These  are  the  hj) 
centrally-ignited,  constant-velocity  rlame 
and  the  Taylor  piston  problem. (2/  For  a 
constant-velocity  flame  Eq.  (3)  reduces  to 


V 


v  ' 3/3  V 

1  .  _JL  Ji 

vbi  vb  Su 
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Therefore,  Fq.  (7)  reduces  to 


p-  =■ 


V 


2y  M1  • 

*  p  [-  r 


i)- 


(8) 


...  where  the  range  on  r  is  Spt  <  r  <  a0t. 

|V'.  -Ecu  at  “ion  (8)  is  identical  to  the  Taylor 
fj  -  solution  for  a  spherical  piston  as  MD  -O' 
(Eq.,  (13) ,  Strehlow  etal.U) ) . 


We  no.te  that  in  these  equations  the,  raft  - 
"of  flame1  ^srea  growth  is  prooortional  to  r 
or-  the  rate  of  decay  of  a  spherical  acoustic: 
Dulse  is,  proportional  to  r-1.-  Thus, ..for  the 
sonercal  case, /the  acoustic  pressure  remains  .., 
'const&nt  at  the  flame  front  and  the  wave 
triaf  is  generated  during,  flame  motion  is  self¬ 
similar  in  r/t.  This  is  a  fundamental  acous¬ 
tic  property  of  a  constant-velocity  flame 
that  is  completely  surrounded  by  combustible 
material  when  it  burns.  It  has  been  verified 
by  numerous  theoretical  and  numerical  invest¬ 
igations.  ,/  . 


APPLICATION  TO  FREE.  CLOUDS 


•Je  restrict  ourselves  to  clouds  of  yf  " 
arbitrary  shape  in  free  space.  To  sinulify 
the  treatment,  we  will  consider  free  clouds 
which  have  a  horizontal  .plane  of  symmetry 
for  cloud  shape,  cloud  distortion  during 
deflagration  and  deflagration  wave  shape. 

Since  we  now  wish  to  simply  examine  the 
implications  of  Eq.  (3)for  general  geometries, 
we.  for  the  moment  restrict  ourselves  to  con¬ 
sidering  only  constant  velocity  flames. 

This  means  that  Eq.  (3)  reduces  Eq.  (4), where 


Su  is  an  effective  normal  burning  velocity 


and  r  is  the  distance  of  an  observer  from 
, . .  the  cloud,  center (  . 


In  general ,  the  evaluation  of  the 
dAf/dt  term  in  Eq..  (4)  -is  complicated 
because  during  deflagrati ve  combustion  the 
cloud  ahead  of  the  flame  will  move  and  ■ 
distort  with  time.  This  distortion  will 
cause  Af  tq  be  a  different  function  of  time 
then  that  calculated  from  purely  geometric 
considerations  using  a  constant  velocity 
flame  propagating  through  a  stationary 
Cloud, 


In  order  to  determine  how  to  best 
aporoximate  the  effect  of  such  distortion 
tor  arbi trary  ignition  of  an  arbitrary  cloud 
one  must  *irst  examine  flow  behavior  in  two 
fimole  1  imit  geometries;  strictly  one- 
dinensional  t low  caused  by  deflagration  and 


volumetric,  heat  ..addi  tion  ‘  in  three  dimens  ions . 


The  wave  diagram  for  a  strictly  one-0 
flow  driven  by  a  conscant-veloc'i by  de'Iacra- 
tion  wave  is  shown  in  Figure  1.  If  we  mooel 
the  combustion-  o‘‘  &  stoichiometric  methane-  :  ' 
si”  mixture  as  proposed  jn  Strehlow  et  al 
the  gas  to  the; rf gVit  of  the  flame  and-  to  the 
lpft  of  the  contact  surface  (regions  0,  1 , 

&  3)  tan.  be  assumed  ,  to  have  fa/' beat  capacity:  ... 
ratio  Yg  -  Y:i  -  y3  =  t. ‘4, ..while  the  gas  jir; 
region  2  has  1.202..  Additionally, 

the  effective  value  by  dimensionless- energy 
addition  is  q  -  8.430.  This  means  (from, 
Appendix  A)  that  6  =  9.242  and  that  Vb/Vu  = 
7.602  for  constant  pressure  combustion. 
Furthermore,  for  this  system  the  Chapman- 
Jouguet  deflagration  Mach  number  is  0.16631. 


'The  properties  of  the. wave  system 
of  Figure  1  were  calculated  using  the  shock 
relationships  from  Liepman  and  Roshko^6) 
for  the  shocks  ST  and  S3,  the  deflagration 
relationship  from  Strehlow  et  al.(3)  fnr 
the  fl ame  and  the  contact  surface  rtqui re- 
ment  that  Up£  =  Up  and  p£  s  P3<  jbe 

results  of  this  calculation  that  are  impor¬ 
tant  to  this  paper  are  summarized  in 
Figure  2. 


Figure  2  shows  that  particle  velocities 
generated  by  flame  propagation  are  relatively 
high  and  that  the  apparent  flame  speed 
relative  to  an  observer  is  from  4.3  to  9.3 
times  the  actual  normal  burning  velocity. 

Also  it  is  less  than  the  space  velocity  of 
a  spherical  flame  if  the  one-dimensional 
flame  is  traveling  at  less  than  about 
0.9  of  its  maximum  (CJ)  velocity.  Further¬ 
more ,  the  ratio  | Up  /Up  !  lies  between  1.0 
and  1.4  in  this  same  velocity  range.  This 
last  observation  simply  means  that  even 
relatively  high  speed  flame  propagation 
from  a  free  surface  in  one-dimensional 
space  displaces  the  gas  rather  uniformly 
to  the  front  and  back.  This  is  to  be 
constrasted  to  the  behavior  of  a  one- 
dimensional  CJ  detonation,  which  displaces 
the  gas  primarily  in  the  direction  of 
detonation  propagation. 


We  now  wish  to  determine  what  happens 
if  we  relax  the  confinement  normal  to  the 
flame  propagation  direction  for  the  rase 
of  deflagretive  combustion.  To  estimate 
this  effect  we  assume  that  if  energy  is 
added  relatively  slowly  to  a  gaseous  volume, 
the  volume  will  expand  an  equal  di s tance 
in  all  directions.  During  'lame  propagation 
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sus^iKiiBsat^^  x  v^«*‘-Vfc'. 
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through  a  cloud  the  flame  can  be  assumed  to 
be  a  relatively  thin  sheet  of  energy  addition 
when  compared  to  the  cloud  dimensions.  In 
other  words  the  flame  can  be  assumed  to 
occupy  a  volume  Aal  when  Al  <<  vR.  Under 
these  circumstances  one  would  expect  the 
major  expansion  to  occur  in  a  direction 
normal  to  the  local  orientation  of  the  flame 
even  though  the  pressure  is  propagating  almost 
equally  in  all  three  dimensions. 

This  assumption  has  been  qualitatively 
verified  experimentally  by  observing  flame 
propagation  in  a  propane-air  elongated  pan¬ 
cake  cloud  above  a  flat  plate  of  about  50  mm 
height  and  0.6  m  x  0.3  m,  edge  ignited  at 
the  center  of  the  0.3  m  edge.^  The  flame's 
apparent  propagation  speed  was  about  3.5 
times  the  normal  burning  velocity,  and  its 
luminous  height  was  about  equal  to  the  initial 
cloud  heignt.  Furthermore  the  flame  flashed 
to  encompass  an  area  considerably  larger 
than  the  original  lateral  extent  of  the  cloud.' 
The  assumption  has  also  been  reasonably  well 
verified  by  numerical  calculations  made  by 
Dr.  Len  Hazel  man  (§)  at  Lawrence  Liverniore 
Laboratories  using  a  2  0  hydrocode.  He 
calculated  flow  and  flame  propagation  in  an 
open  2  D  channel.  His  calculations  yielded 
a  flame  height  almost  twice  the  initial  cloud 
height  and  an  apparent  flame  speed  about  three 
times  the  normal  burning  velocity.  However, 
his  numerical  flame  had  a  thickness  only  one- 
third  of  the  cloud  height  and  this  should 
lead  to  more  vertical  distortion. 

These  results  indicate  that  for  any 
arbitrary  free  cloud  and  igniter  geometry  in 
which  the  flame  is  not  completely  surrounded 
by  a  combustible  mixture  one  nay  make  the 
following  conservative  assumptions  about 
fl ame- induced  distortion. 

1.  The  burning  velocity  that  should  be 
used  is  Su.  This  is  because,  irrespective  of 
the  motion  of  the  flame  induced  by  the  flow, 
the  entire  cloud  will  burn  in  a  time,  t  , 
given  by  Sut  =  rm  where  rm  is  the  distance 
from  the  igniter  to  the  farthest  point  in 
the  initially  quiescent  cloud. 

2.  The  increase  of  flame  area  in  a  free 
cloud  due  to  distortion  of  the  cloud  by  flame- 
induced  flow  will  never  exceed  the  value 

V b/^u  calculated  for  isoboric  combustion. 

This  is  because  this  is  the  value  of  distor¬ 
tion  of  the  unburned  material  produced  during 
combustion  of  a  completely  surrounded  (sDher- 
ical)  flame,  and  relief  due  to  product  gas 
motion  away  from  the  flame  can  only  reduce 


this  factor.  . 

Thus  the  general  procedure  for  determin¬ 
ing  a  first  conservative  approximation  to  the 
effect  of  cloud  shape’ and  .igniter,  location 
on  the  overpressure-time  henavior  in  the 
blast  wave  when  the  igniter  is  not,  completely 
surrounded  by  the  cloud  is  as  fellows: 


1. 

?  . 


4. 


Assume  a  value  cf^$u. 

Calculate  dAf/dt:'.  for  a  flame 
traveling  at  velocity  $u  tnrough 
the  quiescent  cloud. 

Correct  tor  distortion  by  multiplying 
V  b/V  u  ‘ 


Examples  for  Specific  Geometries 

a.  Central  (line)  ignition  of  a  pancake 
cloud. 

The  flame  area  is  given  by 
Af  **">«.■ 

where  rf  is  the  flame  radius  and  H  is  the, 
height  of  the  pancake.  Therefore 

dA,  dr., 

dt  ■  dt  u 


and 


as3- 


max  "  2afR 


w 

tri 


O  ,OOS  l  UJ 


(9) 


where  Robs  is  the.  .distance  of  an  .observer 
from, the  center  of  the  cloud. 


b.  Edge  (line,'  ignition  of  a  oancake 
cloud. 


The  area  of  a  flame  orooagating 
from  an  edge  of  a  cloud  of  radius  R  is  given 
in  terms  of  the  distance  propagated  by  the 
formul a 


Ac  ~  2Hrf cos"1 


r . 
r 

2R 


Now 


dA.  dA, 


__r  s 
dt  "  dr*  dt 


dA, 

dr7 


it"  rfiHlLvti.-'.aAii;" 


i.i 


Therefore 


5  * 


COS  •  nfT  - 


r 

2R 


r 

2^ 


-  f 


V  1  - 


-£j 

2Rj 


or 


max 


$s;h 

4aoRobs 


rv. 

■  !r 

l  u. 
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Hnich  is  exactly  half  that  predicted  for  a  , 
.centrally-ignited  pancake. 

c.  End  (point)  ignition  of  an  ellipsoid 
(cigar-shapcd)  cloud. 

Assume  a  flat  flame  of  location  rf 
propagates  from  one  end  of  a  i  or.g  cloud  with 
a  minor . radius  of  rotation  d  and  major  radius 
b. 

Flame  area  is  therefore 


Af  =  md2 


-  p-j 


Vs  ’*  Vp  *  vr,  where  the  subscripts  s,  p,  anc 
c  refer  to  spherical  ,  oar.cake  and  cigar 
shaped  clouds,  respectively.  Also  define 
the  aspect  ratio  of  the  clouc  as  *>.0  - 
2Rp/Hp  and  A>,Q  =  Dc/Lc.  Using  this 
notation  we  obtain 


H° 


and 


D  = 


raj 


02) 


03) 


Combining  Eqs.  (5),  (9),  (10)  and  (11)  with 
Eqs.  (12)  and  (13)  yields 


max 


12 


—  -  l]  /ft 
v 

'o  : 


1/3 


04) 


and 


or 


dAf  =  “  *">  {£ 


Now  let 

“  b  »  -  L/2  <  r  <  1/2  =  b  and  0  *  2d 


where  L  is  the  length  of  the  ellipsoid. 
Therefore 


fVi 


p  = 


2aJH7  rf 


and 


max 
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Now  assume  that  3,  Su,  a0,  and  (Vl/vJ 
are  Independent  of  cloud  shape  and  that  trie 
observer  is  at  a  fixed  distance  from  the 
center  of  a  cloud  of  fixed  volume  but  of  these 
different  shapes.  Call  the  volun.es 


nax 


-  + 1[ 
0  J 

Q 


i/s 


(15) 


where  the  factor  +  1  =  Vu/Vy  arises 
because  we  are  considering  tne  initial 
cloud  diameter  and  Eq.  (5)  is  written  in 
terms  of  the  final  radius.  Eq.  (14)  is  for 
central-line  ignitior^  of.  a  Dancake:  for 
edge-line  ignition  (Pc/ps^max  is  one  ha1f 
of  this  value. 

Figure  3  is  a  olot  of  Eqs.  (13)  and 
(14).  It  shows  that  the  deflagration, 
after  edge  ignition  of  a  large  aspect-ratio 
cloud,  produces  much  lower  overpressures 
than  central  ignition  of  a  spherical  cloud. 
It  should  be  DOinted  out  that  ooint  ignition 
of  the  edge  of  a  free  spherical  cloud 
>RC  *  1  would  yield  a  higher  overpressure 
than  calculated  with  Eq.  (14)  because  the 
derivation  of  Eq.  (11)  assumed  a  flat  flame. 
In  reality  one  would  ^xoect  a  flame  of 
roughly  sohericai  shape  to  propagate  away 
from  a  point  source  and  the  initial  rate 
of  flame  area  increase  would  be  larger  in 
this  case  than  for  a  flat  flame. 
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effect  of  cloud  size 

* 

Equation  (3)  can  be  used  to  give  an  order 
of  magnitude  estimate  of  flame  areas  and 
accelerations  that  are  necessary  for  a  def¬ 
lagration  to  produce  a  damaging  blast  wave. 

If  we  define  the  threshold  of  damage  to  be 
Bar  and  assume  that  this  level  of  over¬ 
pressure’  is  to  be  produced*! 00  m  from  the 
cloud  center  in  an  atmosphere  with  a  velocity 
orTertmcT5'f~350  m/s,  Eq.  (3)  reduces  to 

dSu  dAf 

'•6  x  10  ’  *f  3T  *  su  dT  06) 

Table  I  was  constructed  from  Eq.  (16)  under 
the  two  limit  assumptions  that  (1)  dS  /dt  =  0 
and  (2)  dAf/dt  =  0.  u 

Notice  from  Table  I  that  even  for  very 
high  flame  velocities  the  rate  of  flame  area 
increase  must  also  be  very  high  if  even  a 
weak  blast  wave  is  to  be  generated.  On  the 
other  hand,  only  extremely  large  initial 
flame  areas,  exhibiting  very  large  flame 
acceleration  over  their  entire  frontal  area 
are  necessary  to  produce  a  weak  blast  wave 
from  deflagrative  combustion.  Both  of  these 
observations  imply  that  extremely  large  clouds 
are  required  if  one  Is  to  produce  significant 
overpressure  by  deflagrative  combustion  alone. 

DISCUSSION  ANO  CONCLUSIONS 

The  application  of  simple  source  acoustic 
theory  to  deflagrative  combustion  of  an  un¬ 
confined  cloud  shews  that: 


1)  It  Is  very  difficult  to  produce  a 
damaging  blast  wave  by  aeflagrative 
combustion  once  the  Flame  is  not  comoletely 
surrounded  by  a  combustible  mixture. 

2)  The  maximum  overoressure  produced,  other 
circumstances  being  equal,  is  proportional 
to  the  ratio  of  the  minor  dimension  of 

the  cloud  to  the  distance  to  the  observer. 

3)  Blast  pressure  is  rather  uniformly 
distributed  in  all  directions,  i.e.  ,  the 
blast  is  roughly  spherical. 

4)  The  cloud  must  be  very  large  if  a 
damaging  blast  wave  is  to  be  produced. 

5}  Spherical  flame  propagation  calculations 
such  as  those  of  Kuhl  et  aj..  and  Strehlow 
et  ah  greatly  overestimate  blast  pressures 
from  deflagrative  combustion  following  edge 
ignition  of  clouds  with  large  aspect  ratios. 

From  various  accident  accounts  of  the 
sequence  of  events  that  led  to  the  production 
of  a  damaging  blast  wave  after  delayed 
ignition  of  a  massive  spill  of  combustible 
material  it  appears  that: 

1)  There  is  a  threshold  spill  size  below 
which  blast  damage  does  not  occur.  Gugan's'*' 
documentation  of  incidents  shows  that  blast 
damage  has  been  observed  for  spills  of  less 
than  2000  Kg  but  mere  than  100  Kg  only  for 
the  fuels  Hg,  H?-C0  mixture,  CHj  and  CoH^. 
Blast  damage  has  been  recorded  only  for 
spills  greater  than  2000  Kg  for  all  other 
fuels. 

2)  In  the  majority  of  cases  where  blast 


TABtE  I 


(1)  dSu/dt  =  0 

(2) 

dAf/dt  =  0 

S.  rn/s 

0-1  CL 

Af ,  m2 

dSu 

u 

t 

m2/s 

dt  ’ 

m/s2 

1 

1.7 

X 

107 

100 

1.7  x 

10s 

10 

1.7 

x 

10s 

10,000 

1.7  x 

103 

100 

1.7 

X 

10s 

1  ,000 ,000 
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damage  occurred,  fire  was  present  for  a 
considerable  period  before  the  blast  occurred. 
3)  In  many  cases  damage  is  hiqhly  direction¬ 
al. 

These  observations  when  coupled  with 
the  results  of  simple  source  acoustic  theory 
for  deflagrative  combustion  lead  to  the 
following  conclusions. 

1.  There  should  be  a  size  threshold  below 
which  blast  will  not  occur  as  long  as 
ignition  is  "soft",  i.e.,  does  not  directly 
trigger  detonation. 

2.  The  fact  that  fire  is  present  early  after 
ignition  indicates  that  massive  flame  acceler¬ 
ations  are  necessary  to  lead  to  blast  wave 
formation.  Since  the  flame  must  have  burned 
through  the  cloud  edge  by  the  time  the  blast 
is  produced,  simple  acoustic  source  theory 
must  be  operative  for  even  high  deflagration 
velocities.  Thus,  the  blast  must  arise  from 
some  sort  of  effectively  supersonic  combustion 
process  or  from  very  rapid  increases  in 
effective  surface  area  of  the  flame. 

3.  Simple  acoustic  source  theory  for  detla- 
grative  combustion  shows  that  deflagrative 
orocesses  per  se  cannot  produce  highly 
directional  effects.  However,  it  is  well 
known  that  detonative  combustion  of  a  cloud 
does  produce  highly  directional  blast  wave 
effects. 
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APPENDIX  A 

The  Generalized  Heat  Addition  Hugoniot 

Strehlow  et  aJL  (3.)  have  shown  that 
a  generalized  heat  addition  Hugoniot  quite 
accurately  fits  the  real  Hugoniot  for  hydro¬ 
carbon-air  combustion  over  the  range  of 
interest  for  vapor  cloud  exDlosions.  Briefly 
they  assume  that 


and 

h!  *  V  -  5 

where  Q  is  the  energy  added  at  the  deflagra- 
ti ve/detonative  discontinuity  and  Cp0  and 
Cp^  are  different  heat  capacities  for  the 
reactive  and  product  gases  respectively. 

Using  this  approach  they  define  a 
dimensionless  energy  addition  as 

3  =  (r,  -  1  )Q/P0V0 
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'Thus  the  volume  change  for  constant  oressure 
combustion  is  given  by  the  expression 
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therefore 
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Now  define 
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for  ease  of  nondimensional  izing  the  acoustic 
problem.  This  yields 

*  .  Yo(YrYo}  .+  !o  a 
4  W1}  Yi 

For  stoichiometric  methane-air  mixtures ^ 

Y-|  =  1.202  and  q  =  o.43G 

Therefore 
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for  a  stoichiometric  methane  cloud. 


Figure  i.  The  wave  diagram  for  strictly  one-dimen¬ 
sional,  constant-velocity  deflagration  of  a  free  cloud  after 
edge  ignition  (edge  located  at  x  =  0,  Sf  is  apparent  flame 
speed). 


Figure  2.  Flow  velocities  and  apparent  speeds  associated 
w  1th  '  strictly  one-dimensional  deflagration  combustion 
illustrated  in  Figure  I. 


* 

Figure  3.  Effect  of  the  aspect  ratio,  .41,  on  the  maximum 
blast  wave  pressure  rise  for  the  deflagrative  combustion  of 
pancake  and  cigar-shaped  clouds.  Cloud  volume,  normal 
burning  velocity  and  observer  distance  from  cloud  center 
all  assumed  to  be  constant  from  cloud  to  cloud. 
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